The Newtonian limit for perfect fluids 

Todd A. Oliynyk *t 
Max-Planck-Institut fiir Gravitationsphysik 
Am Muhlenberg 1, D- 14476 Golm, Germany 



Abstract 

We prove that there exists a class of non-stationary solutions to the Einstein-Euler equations which 
have a Newtonian limit. The proof of this result is based on a symmetric hyperbolic formulation of 
the Einstein-Euler equations which contains a singular parameter e = vt/c where vt is a character- 
istic velocity scale associated with the fluid and c is the speed of light. The symmetric hyperbolic 
formulation allows us to derive e independent energy estimates on weighted Sobolev spaces. These 
estimates are the main tool used to analyze the behavior of solutions in the limit e \ 0. 



1 Introduction 

The Einstein-Euler equations or, in other words, the Einstein equations coupled to a simple perfect fluid 
are given by the following system of equations 

G'^ = (1.1) 

V,T*^' = (1.2) 

where the stress-energy tensor for the fluid is given by 

T*^ = {p + c-^pyv^ + pg'^ (1.3) 

with p the fluid density, p the fluid pressure, and v the fluid four-velocity normalized by v^Vi = — c^, c the 
speed of light, and G the Newtonian gravitational constant. The study of the behavior of solutions to 
these equations in the limit that e — vt / c \ where vt is a characteristic velocity scale associated with 
the fluid matter is known as the Newtonian limit. By suitably rescaling the gravitational and matter 
variables (see section [2]), the Einstein-Euler equations can be written as 

Qij ^ 2e4r*i and V,T^^ = (1.4) 

where Viv"^ = —e^^, and t — x'^/vt is a "Newtonian" time coordinate. In the limit e \ 0, one expects 
that there exists a class of solutions to Einstein-Euler equations (|1.4p that approach solutions of the 
Poisson-Euler equations 

dtp + diipw') = 0, (/,J = 1,2,3) (1.5) 

p{dtw^ + w'diw^)^-{pd''<i> + d-^p), {d' = 5'-'dj) (1.6) 

A$ = p, (A = 9/9^) (1.7) 

of Newtonian gravity in some sense. As above, p and p are the fluid density and pressure, respectively, 
while is the fluid (three) velocity. This problem has been studied since the discovery of general 
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relativity by many people and there is a large number of results available in the literature. The majority 
of results are based on formal expansions in the parameter e which are used to calculate the (approximate) 
values of physical quantities and also to investigate the behavior of the gravitational and matter fields in 
the limit e \ 0. For some classic and recent results of this type see [2,3,6,9,13,20 22,31,41] and reference 
cited therein. The main difficulty with the formal expansions is that they leave completely unanswered 
the question of convergence. In the absence of a precise notion of convergence, it becomes unclear to 
what extent the formal expansions actually approximate rclativistic solutions. 

In this paper, we go beyond formal considerations and supply a precise notion of convergence for 
gravitating perfect fluids as e \ 0. This necessitates introducing suitable variables that are compatible 
with the limit e \ 0. The metric gij, which defines the gravitational field, turns out to be singular in 
this limit. To remedy this problem, we introduce a new gravitational density u^^ which is related to the 
metric via the formula 

(1.8) 



V-det(Q) 
where 

From this, it not difficult to see that the density u'-' is equivalent to the metric g^j for e > and is well 
defined at e = 0. For the fluid, we also introduce a new velocity variable according to 

v^=w' and v^ = l + €W^. (1.10) 

For technical reasons, we only consider isentropic fiow where the pressure is related to the density by an 
equation of state of the form p = f{p). Moreover, to formulate a symmetric hyperbolic system for the 
fluid variables {p, f}, we need to deal with the well known problem that the system becomes singular 
when p + c~^p = 0. This is a particular problem for fluid balls having compact support. To get around 
this problem, we follow Rendall [34] and use a technique of Makino [24] to regularize the fluid equations 
so that a class of gravitating fluid ball solutions can be constructed. Thus as in [34], we assume an 
equation of state of the form 

J5 = (1.11) 
where K e ]R>o, n gN, and we introduce a new "density" variable a via the formula 

As discussed by Rendall, the type of fluid solutions obtained by this method have freely falling boundaries 
and hence do not include static stars of finite radius and so this method is far from ideal. However, in 
trying to understand the Newtonian limit and post-Newtonian approximations these solutions are almost 
certainly general enough to obtain a comprehensive understanding of the mathematical issues involved 
in the Newtonian limit and post- Newtonian approximations. We would also like to remark that the 
results contained in this article are largely independent of the specific structure of the fluid equations. 
We therefore expect that the analysis in this paper can be carried over without much difficulty to other 
matter models whose equations can be formulated as a symmetric hyperbolic system and have a finite 
propagation speed for the matter density in the limit e \ 0. 

Our approach to analyze the limit e \ is to use the gravitational and matter variables {u'-' , w%a} 
along with a harmonic gauge to put the Einstein-Euler equations into the following form 

b°{eV)dtV = ^c'diV + b\e, V)diV + f{e, V)V + ^g{V)V + h{e) , (1.13) 

where V comprises both the gravitational and matter variables, and the are constant matrices. This 
system is symmetric hyperbolic and hence by standard theory there exists local solutions. However, 
the difficulty in analyzing the limit e \ of such solutions is that the equation contains the singular 
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terms e~^c^diV and e~^g{V)V . Although, singular limits of symmetric hyperbolic equations have been 
previously analyzed in [5,19,37,38], these results cannot be directly apphed to the system (|1.13|) . There 
are two main difficulties in adapting these results to the Einstein-Euler system. The first is that the 
Einstein-Euler system (j6.ip must be modified by including an elliptic equation, essentially the Newtonian 
Poisson equation, in order to be of the canonical form required by [5, 19,37,38]. This results in a coupled 
elliptic-hyperbolic system of the form 

B°{eW)dtW = ^cJdiW + S^(e, W)diW + F(e, W)W + H{e) , (1.14) 

where W is related to V via an elliptic equation and F is a non-local functional. The second difficulty is 
that the initial data which must include a l/r piece for the metric and hence it cannot lie in the Sobolev 
space . This l/r type fall-off behavior is crucial for obtaining the correct limit and is intimately tied 
to the elliptic part of our formulation of the Einstcin-Eulcr system. The standard procedure in general 
relativity to deal with this type of fall off, at least for elliptic systems, is to replace the spaces with 
the weighted Sobolev spaces [1,7]. However, the arguments used in [5, 19,37,38] fail for the weighted 
spaces as the weight used to define the Hg spaces destroys the integration by parts argument which is 
used to control the singular term e^^c'djW in (|1.14p . Indeed, using integration by parts, it follows easily 
from the definition of the weighted Lg inner-product (see (jA.4p with e = 1) that 

(~e-'c'diW\W)L2 = -l(djia-^'-^)c'W\W)L2 (1.15) 
2e 

where a{x) ~ \/l + \x\'^ /A. In general, this term will in blow up as e \ unless 6 ~ —3/2 which coincides 
with the standard norm. However, to include l/r fall-off, we need to consider —1<S<0 which 
introduces a singular 1/e term into energy estimates based on the weighted norm Hg. 

To overcome this problem, we introduce a sequence of weighted spaces Hg^ (see appendix \X\ for a 
definition) by replacing the weight (t{x) with (Tg{x) — a{ex). Under this replacement, ()1.15p changes to 

{-e-'c'diW\W)Li^ < C{W\W)l2^^ , 

which is no longer singular as e \ 0. This allows us to derive e independent energy estimates for 
solutions to the Einstein-Euler equations. These estimates can then be used to define a precise notion 
of convergence for gravitating perfect fluids solutions in the limit e \ which is essentially a statement 
about the validity of the zeroth order expansion in e. This is formalized in the following theorem; for a 
more precise version see propositions 15.11 16.11 and 17.81 and theorems 17.71 and 17.121 

Theorem 1.1. Suppose -1 < 5 < -1/2, k >3 + s, jS^ £ flLo C'^dO' ^1' -^^5-1 ) « harmonic gauge 
source function, and a,w^ G Hg_-^^, }^'^ G H^^^ , 34'^ G Hg_^ is the free initial data for the Einstein-Euler 
equations where supp a C B*„ for some R* > 0. Then for small enough, there exists a T G (0, T*] 

o 

independent of e £ (0,eo], and maps 

s+l 

u:^(i)-u^(0), div^^it), dMt), a,(i), wlit) ef]c\[0,nHtU 

$ G C°{[0,T*],H^+^) n C\[0,T*],H^+^) , 
u;^ G C° ( [0, T*] , Ht 1 ) n ( [0, T*] , H^l ) , 
p G C0([0, T*],Hl,) n C\[0, T*],Hlzl) , 

such that 
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o 



ag(0) — a , 

o 

p(0) = p = {4Kn{n + l))-"a^" , 

o o 

where (j)^ = 4>{e, PjW^ (3^ {0),^^"^), and = ^{e, p,w^ , f3^ {Q),^^"') is the initial data deter- 

o o o o 

mined by the gravitational constraint equations ( see proposition \5.1\) . and gij is determined from 

\x\i{0) by the formulas (dH]) and ([XT]) . 

(a) {u^J {x^ ,t),ae{x^ ,t),wl{x^ ,t)} determines, via the formulas (jl.Sp . (|1.9p . (jl.lOp . anrf (|1.12p . a 1- 
parameter family (0 < e < eo) o/ solutions to the Einstein- Euler equations (jl.4p m f/ie harmonic 
gauge edtuf^ + djvJJ — ej3^ on the common spacetime region {x^ , t) E D = M.^ x [0, T], 

(Hi) {^{x^ ,t), p{x^ ,t),w^ {x^ ,t)} solves the Euler-Poisson equations (jl.5p - (|1.7p on the spacetime region 
D, 

(iv) there exists a constant R £ (i?*,oo) independent of e £ (0,eo] such that supp a^it), supp p{t) C 
Br for all (i,e) G [0,r] x (0,eo], and 

(v) there exists a constant C > independent of e E (0,eo] such that 

iiu:^(i) - <5^<s^$(i)iiL« + II - 5\5idMmH>^-^ + \W{t) - w'{t)\\H.-i 

+r^\\v^{t) - 1\\h^-i + \\p,{t) - p{t)\\H^-i + \\dtp,{t) - dtp{t)\\H^-2 < Ce 
for all {t,e) G [0,T*] x (0,eo]. 

We remark, that the techniques of this paper can also be used to derive convergent expansions in e 
of the type considered in theorems 2 and 3 of [19] and [38], respectively. These convergent expansions 
in general differ from the formal post-Newtonian expansions. To get post-Newtonian expansion to a 
certain order in e requires that the initial data must be chosen correctly. In the absence of constraints 
on the initial data, a general procedure for doing this is discussed in [5]. Due to the fact that there 
are constraints on the initial data in general relativity, this becomes a non trivial problem called the 
initialization problem. See [18] for an extended discussion. The proof of convergence and a discussion of 
the initialization problem will be presented in a separate paper [27]. 

We note that similar results for the Vlasov-Einstein system have been derived in [36] using a zero 
shift maximal slicing gauge. However, unlike [36], our approach is able to handle not only higher order 
expansions in e, but also a wide variety of matter models. We also note that in [16, 18], there is another 
interesting proposal for analyzing the limit as e \ which is based on a gauge for which the Einstein 
equations are again elliptic-hyperbolic but distinct from [36]. As in this article, the authors of [16,18] 
also propose to use the methods of [5, 19,37,38]. However, the required estimates are not proven and it 
is yet to be verified if this approach would be successful. 

We remark that the results of this and the companion paper [27] are local in time and therefore 
address the "near zone" problem. In the special case of spherical symmetry, the situation improves and 
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there are some global results available on the Newtonian limit [26,32]. However, because spherically 
symmetric systems do not generate gravitational radiation, these results do not shed light on the "far 
zone" problem for post-Newtonian expansions where radiation plays a crucial role and the e \ limit 
must be analyzed in the region "close" to future null infinity. We plan to investigate the far zone problem 
in the near future. 

Our paper is organized as follows: in section [51 we define dimensionless variables for the Einstein- 
Euler system. Sections |3] and |4] are devoted to introducing variables and a gauge condition that cast 
the Einstein-Euler equations into a form suitable for analyzing the limit e \ 0. Appropriate initial 
data which is regular in the limit e \ is constructed in section [S] while in section [5] we prove a local 
existence theorem for the Einstein-Euler system on the weighted spaces. Finally, in section [71 we show 
that solutions to the Einstein-Euler system converge as e \ to solutions of the Poisson-Euler system of 
Newtonian gravity. A precise statement of convergence is contained in theorem 17.121 which is the main 
result of this paper. 



2 Units 

Our conventions for units are as follows: 

W]=L, = [p]-^, H = [v']^[c] = ^. and [G] = 



Note that with these choices the stress-energy tensor has units of an energy density, i.e. [T^^] = To 
introduce dimensionless variables, we define 

= vtv^ and p — pxp 

where vt and px are "typical" values for the velocity and the density, respectively. The Einstein-Euler 
equations then can be written as 

G''J = 2e^f'^ and = 

where 

_Vt _ AnGpT ~i _ /- i . _ - _ P 

€ — , K — „ , X — \/ K X , gij — gij , p — „ , 
c vj. v^pt 

and 

f'^ = ip + e^pyv' +P9"^ ■ 
The normalization u^w* = — c^, implies that 

WiW* := gijv'v^ = . 
Also, we can introduce a time coordinate t via 

t = X^/VT ■ 

With these choices, we have 

[e]^[v'] = [p] = [p] = [g]^[x^]^l, M = ^, [t] = [T], and M - ;^ • 

Thus all our dynamical variables and coordinates are dimensionless and the two constants vt and k can 
be used to fix the length and time scales by using units so that 

Vt — i and k ~ 1 . 

In this case we can use t and x^ interchangeably as long as we remember that they carry different units. 
To simplify notation, we will drop the "hats" from the hatted variables for the remainder of this article. 
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3 Reduced Einstein Equations 

To aid in deriving the appropriate symmetric hyperbolic system for the gravitational variables, we tem- 
porarily introduce a new set of coordinates related to old ones by the simple rescaling 

x'^ = x'^ , x'^^x^/e 

and let 

9. = — a. = — 

dx^ ' dx^ 

In the new coordinates, the metric g^j and its inverse g'^ are given by 

Next, consider the metric density 

r = /^3*^' where |g| = - det(5.,) . (3.2) 
We note that the metric g"^^ is related to the density by the following formula 

g'^ = -^q'^ where I5I = -detg*^' , (3.3) 

and hence 

ia'^) = ^ilL (3-4) 



To obtain a gravitational variable that is regular and non-trivial in the limit e \ 0, we define 

^=^^2(0''-^'') (3-5) 

where 



-1 



is the Minkowski metric density. As stated in the introduction, for e > 0, the metric gij can be recovered 
from the density ii*-' via the formulas (|1.8p - (|1.9p . As we shall see, even though the metric gij is singular in 
the limit e \ 0, the quantity \x^^ is well defined at e = 0. We note that these variables are closely related 
to the gravitational variables discovered by Jiirgen Ehlers and subsequently used in the papers [17,28,29] 
to construct stationary/static solutions to the Einstein equations coupled to various matter sources. 
In the (x*) coordinate system, the Christofell symbols are given by 

f% = £2 ^-fc™(20,,g^p - -Q,fQ,p)d^v.'P + 2{Q,p5\^d,)u'P - 20,(,9,)u")) . (3.6) 

We note that Christofell symbols in the (a;*) coordinate system are related to the T^j as follows 

^44 = £ 44 , r44 = e I4 , r^4 = f j^4 , (3.7) 
T\g — eT\g , r^4 = e ;^4 and Tg^^ = T^c ■ (3-8) 

Using (|3.6[) . a straightforward calculation shows that the Einstein tensor G*-' is given in terms of the 
density u*-' by 

gv ■.= ^\g\ G'^ = fl^^fcfS'^' + [A'i + B'^ + C'^) + D'^ (3.9) 
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where 

|5| =-det(r"), (3.10) 

A'' = 2(iflfe,0„„ - 0fc„0£„) (fl'Pfl^^ - ifl*^flP9)apu"a,u'"" , (3.11) 

B'^ = 40M(2r''^™u^)^a„u'=™ - ir^5„u'="5„S"^ - 0'""a™u*'=a„uJ"^) , (3.12) 

C'^ = 4(afea*^5£u" - Bku'^Bi^u^^) , (3.13) 

D'^ ~ fl'^S^fu" - 29^^u'=('0^)^ . (3.14) 

To fix the gauge, we assume that 

d,u'^ = efi^ (3.15) 
for prescribed spacetime functions (3^ — (3^{x^ ^x'^). For e > 0, di\x^^ — imphcs that 

d,-0'^ = 4e3/3-'" 

or equivalently 

dkd''^ = Ae^P^ and dks''^ = Ae^f]"^ 



where g*-' — ^ — det{gki)g^-' is the metric density in the {x'') coordinates. Thus (I3.15P is, for e > 0, a 
generahzed harmonic type gauge and is harmonic if the functions (3^ are chosen to be identicaUy zero. 
Clearly, if we define 

E'^ := Q^^dkp" - 2dkpH'^'' , 

then (|3.15p implies that 
Setting 

G'^ G'^ - D'^ + eE'^ = Q^'^Bl^u'^ + eE'^ + e^{A'^ + B'^ + C'^) (3.16) 

and 



the Einstein equations G'^ = 2e^r*-' in the gauge p.lSp become 

G'i = T'' (3.17) 

which we will refer to as the reduced Einstein equations. 

To write the reduced Einstein equations in first order form, we introduce the variables 



diu'^ if fc = / 
eB^u'^ if fc = 4 



The reduced Einstein equations then become 

-Q'^'^Biu'^ = 2g*'Biuii + 3"Biuj + eE'^ + t^{A'^ + B'^ + C'^) - T'^ , 
9"B,iej = 9"Bju^ , 

B4U'' = u:i , 

or equivalently 

e e ^ e 

C4UJ — -0 CJU4 , 



1 

e '4 
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eu^^ <-':=u^^ (3.18) 



Next, define 

■ 

and let 

V = { (r*^) e M4x4| det(?7'^' + 4r'^') > } . 

Then using vector notation 

we can write the reduced Einstein equations as 

A'^{eu)diU'^ = jC'diu'^ + A^{u)diu'^ + F'^ (e, u) - i {T'^ , 0, 0)^ (3.19) 



where 



and 



/I - 4eu*^ 0\ 

A'*(eu) = 5^-^ + 4eu" , (3.20) 

V y 

( 5" 0\ 

C' =\6" , (3.21) 

\ 0/ 

^8u'*^ 4u^'^ 0\ 

^^(u) = I AvJ'' , (3.22) 

0/ 



F'\e,n) = {E'^ + (eu,Ufc),0,<^)^ . (3.23) 



The functions /*-' (eu, Ufc) are analytic for eu e V and moreover are quadratic in Uk- Here we are using 
the notation 

u=(u'^) and Ufe = (4^). 
The stress-energy tensor is given in terms of the u variable by 



and hence 
where 



^^^'^ = P [ig^v^ e-i [{\g\ - l){vr + {{vy - 1)] 



+ '1^1 V ept;^t;4+4e|5|-VV4j ^(t;^)^ + |^|-i/2p(-l + 4eu^4)J • (•^•^b) 

We remark that if — 1 = 0(e), then S^^ is regular in e as is easily seen from the above formula and the 
expansion 

|5| =l + 4er;,,u^^ +/(eu) (3.27) 
where /(eu) is analytic for eu e V and also satisfies f{y) = 0(|yp) as y ^ 0. 
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4 Regularized Euler equations 

There are various approaches to symmetric hyperbohc formulations of the relativistic Euler equations 

[4, 14, 15, 34, 40]. We use the approach of [4] which is based on fluid projection and the introduction of a 
Makino variable. 

In the coordinates (a;*), the Euler equations are given by 

V.T'J = (4.1) 

where T'-' = (p + e^p)v^v^ + pg''^ and the fluid velocity is normalized according to 

v^v' = -\. (4.2) 

Differentiating (|4.2p yields 

vS/jV^ = (4.3) 

which implies 

v'viVjv' = . (4.4) 

Writing out (|4.ip explicitly, we have 

{d,p + e^dipyv' + {p + €^p){v^Viv' + v'ViV^) + g'^d.p = . (4.5) 

The operator 

H = Si + e'^v^v^ 

projects into subspace orthogonal to the fluid velocity v^, i.e. L^Ll. = Ll. and L^v^ = 0. Using L^. to 
project the Euler equations (|4.5p into components parallel and orthogonal to v'^ yields, after using the 
relations (|4.2|) - (|4.4p . the following system 

v'^^p +ip + e^p)L]V^v' = , (4.6) 

M.^v'^VkV^ + -^—VAp = , (4.7) 
p + e^p ■' 



where 

As discussed in the introduction, we introduce a new density variable a via the formula (jl.l2p . 
Multiplying ()4.6p by the square of the function 



gives 



h^v'd^a + h\p + e^p)^L)V,v^ = 0, (4.8) 
dp ■' 



dp 

p + e^p da 

where 



M,,v^^kv' + —-^-fLld.a = , (4.9) 



2 dp 1 2 

■5 = 3~ = T^" 
dp 4n^ 

is the square of the speed of sound. A simple calculation shows that 

s2 dp 2/ , 2 

h [p + e p)— = q 



p + e^p da dp 
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where ^ 

q = q{e,a) = — — — -a. 

2nn(ea) 

This shows that the system (j4.8|) - (l4.9p is symmetric, and moreover at a point where a = and hence 
p = p = 0, it is regular unhke (|4.6p - (|4.7p . This is the point of introducing the Makino variable a. Also 
note that the pressure is given in terms of the Makino variable by 



Define 
so that 

Using vector notation 

we can write (|4.8p and (|4.9p as 

where 



^ (4is:n(n + l))"+i 

I -I J 4 -4 

w := V , and w :— v 

e 

/ „„/ „„J „,4 _ 1 , „,.4 



?; = w 



and = 1 + ew'* . (4.11) 



w = (a, w*)"^, 



a''a4W = a^diw + b (4.12) 



/i2(l + ew4) eqL^ 
eqLf My(l + ew4) 



(4.13) 



and 



(4.15) 



From (031), (jXTH)) . and ([XTT]) . we find that 

ffy =?7»j (4.16) 
where the fij{y) are analytic and satisfy fij{y) — 0(|y|) as y — > 0. Also, (|3.6p shows that 

f^- =e[r7'="(2,7,,,7,p-,yy,7,,)6U^ + 2(,7,p<5f,eu^^-2%(,e4^^ + 6/,^-(eu, eu,„) (4.17) 
for functions /*j (eu, eUm) that are analytic for eu G V, linear in the eUm, and satisfy (0, y) = 0. So then 

My = g,j + 2e^gikgjiv^v'^ = 5ij + m^- (eu, eu;'') (4.18) 

and 

Ll = (5^ + e'g,kv''v' = (5^' - 5^51 + ^^eu, eu.'=) (4.19) 

where ^■^(eu, ew'') and (eu, eiu'') are analytic for eu € V and £^(0,0) — my (0, 0) — 0. Using (|4.16p - 
(|4.19p . the matrices a' and the vector b can be written as 

a^=(^J ^"^.^ +«^eu,ew), (4.20) 
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and 

, / \ / abi{eu, ew) • euA 

V-?7™(2?74£?74p + mpyP - 2(77,p51uf - 277Muf ); ^ ^ ^2(eu, ew) • / ' ^ ' ^ 

Note that (i) a'^, a, , bi, and 62 are analytic in all their variables provided that eu G V, (ii) a^, a 
and are symmetric, and (iii) a*(0,0) = 0, a^(0,0) = 0, a(0,0) = 0, 6i(0,0) = 0, and 62(0,0) = 0. 
Consequently the system (|4.12p is symmetric hyperbolic on a region where (eu, ew) is small enough to 
ensure that is positive definite. This can always be arranged by taking e small enough and since we 
are interested in the limit e \ no generality is lost in assuming this. 

It is important to realize that the derivation above of (|4.12p required that both the Euler equations 
(|4.ip and the fluid velocity normalization (|4.2p are satisfied. Alternatively, we can first assume that (|4.12p 
is satisfied and then show that (|4.ip and (|4.2p are also satisfied. To see this, define 

TV := evy + 1/e = eg^iil/e + w^ f + 1/e + 2,g4j(l + ew^)w' + egijw'w' . (4.23) 

Clearly, J\f — is equivalent to v'^Vi = — l/e'^ for e > 0. Furthermore, any solution of (|4.12p also solves 
(|4.6p - (l4.7p for any e > 0. So assuming that u is a solution to the system (|4.6I) - (|4.7I) . contracting (|4.7p 
with yields 

,2-1- \ -,k i 



For (2e7V - 1) 7^ 0, this imphes 



(1 + ew^)diAf - -w'diAf + ^^^AA . (4.24) 

Clearly, this is a symmetric hyperbolic equation for Af whenever < 1/C < (1 + ew*) < C for some 
constant C. This can always be arranged at x*^ = by choosing e small enough. Therefore, if initially 
A/'|_^4_Q= 0, then A/" = for as long as (1 + ew'^) stays absolutely bounded and bounded away from zero. 
Consequently, choosing initial data for the system (|4.12p such that ■^\^i_Q= will guarantee that the 
solution will satisfy the full Euler equations (|4.5p in an open neighborhood of the hypersurface x'^ = 0. 
In particular, if {a, w'} is a solution to (|4.12p with initial data satisfying ^1x^=0^ then a is a solution to 
the equation 

940! + X^dia + Ya = (4.25) 

where ^ 

:= ^ , , and Y := — , . (4.26) 

Observe that 

Y — Y{ew*, ea){edtw'^ + djuu^) + Y{eu, ew^, eufc, ew^ , e) 

where F(0, 0) — l/(2n) = 0, y(0, . . . , 0) = and F(ew*, ea.^)^ Y{eu, ew*^, eUfc, ew^, ea^) are analytic on the 
region eu e V and 1 + ew'* > 0. 

5 Newtonian initial data 

Let Sq = R'^ be the hypersurface defined by 5*0 := {(a;^, 0) | [x^] S R^}. The covector rii — Sf is conormal 
to 5*0 implying that constraint equations for the initial data on So are given by riiG^^ — 2e*riiT'-' . Defining 

:-e-i(g4J_^4j) C^:=g^i-T^\ 
we find that = is equivalent to UiG''^ ~ 2e'^n.iT'^^ for e > 0. Also, by defining 

W := d,u'' - ef3^ , (5.1) 
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the generalized harmonic gauge (|3.15p can be written as W = 0. 

As will be seen in the proof of the next proposition the equations = are regular at e = 0. So to find 
appropriate initial data that is well defined at e = 0, we solve the regularized constraint equations = 0. 
Moreover, we must also ensure that the harmonic gauge condition W = and the fluid normalization 
Af = are satisfied. To solve the constraints = 0, W = 0, and A/" = 0, we use a implicit function 
technique based on the work of Lottermoser [23] . We assume that the fluid velocity can be written as 
(|4.10p which is consistent with the expected behavior of the fluid velocity as e \ 0. We will not assume 
that the density and pressure are related by the equation of state (|l.lip . Instead, we will consider them 
as independent prescribed fields for the purpose of finding solutions to the constraint equations. We do 
this so that the following proposition remains valid for other equations of state. 

Proposition 5.1. Suppose -1 < (5 < 0, fc > 3/2 + 1, i? > and {p,p,w^ ,l" , f3^ £ [H'^zIY x 
Hs_i X {H^zlY X Bji{Hg) . Then there exists an eo > 0, an open neighborhood U of {p,p,w^ , i^-' ,1 )> 
and analytic maps (-eo,eo) x [/ — > Hg_^ : {e, p,p,w^ PZi^''')^ uj^, (-eo,eo) x U : 
{e,p,p,w' ,(3Zi")^ 4>, (-eo,eo) x U H^^ : {e, p,p,w' , P^i") ^ such that for each 
{PtP^w^ , j3Zl^'^) € U, {e, p,p,w^ jW'^jU^ , l3Zd4U^-') is a solution to the three constraints 



wher 



^0, W ^0, and U ^0, (5.2) 



id^^-n = [ _^Jh ^ ,j I^^CtS) {t = x')^ (5.4) 



-dK^-' + 13'' -dxn,'^ + f3' 



4 1 , -e.94,/w"' - ^e^igijw-'Y - gaie^gijw'w'' + 1) 
w — 1 . (5.5) 

e £344 

Moreover, if we let 0o = 0|e=O; — ^'\c=q, and = w*|e=Oj then 0o, ItIq, and Wq satisfy the equations 
A(j)o = p , Aw'o - da" + pw' , and ^ , 

respectively. 

Proof Let u^'* = 0, u^^^ = ea", u^^ ^ elr^ and ^411^"' = e]^'/ . Solving ^^1^^^= yields 

diU-^^ = e(-9/lr^ + /J^) and Biu'^' = e{-dii" + P'^) (5-6) 
while solving A/"!;, = gives 

I On 



4 1 , -^dijw' - V e-^igijwJ)^ - gaie'^gijw'wJ + 1) 

w = 1 z • (5.7) 

e eff44 

From (|3.3|) and ([33]), it is not difhcult to verify that 

= e-^f{ew^,e^i,e^ni,e^(j3) 

where /(y) (y = (?/i, . . . , yi)) is analytic in a neighborhood of (0, 0, 0, 0) and and moreover /(y) = 0(|yp) 
as y ^ 0. 

Using the relation (|5.6p to eliminate 9411'*^ and 9411'*'^ in favour of tti^ and j^'^, we find that 
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where 

Using this and equation ([XTUl) - ([XTi|) . ([XMjl - f^:^ . and (|iTI]l) - (|ilTl) . we see that 

= Atr^ + dLi" + eh' 

+ efie^i, e^ttj, ^<P, eD^, eDfo, D<l>, £34, ^(-^Ktr^ + f3^), e(~dK}'''' + /3^)) - S^' , (5.8) 

and 

C* = Acb-p- e{2w^ + e(u>4)2)p - e^i.^^"^^ + 4e2/i4+ 

e'tt', e^tr, 534, e(-9A'ft)^ + /3^), £(-9^3^^ + - e^^^ , (5.9) 

where the functions /■'(y) (y = (j/i, • ■ • , J/g)) are analytic in a neighborhood of {(0, 0, 0)} x U where U is 
any open set and are quadratic in (2/4, . . . , j/g). Note that 

= pStit, w', e^i, e^Xo, e^) + p5f (e, u.^ 6^3, e^xo, e^) 

and 

= pw' + epS'((e, w^, e^3, e^tr, e^) + epS'|(e, w^, £^3, e^ro, e0) 

where the functions S^{y) (y = (yi, . . . ,2/7)) are analytic in a neighborhood of [/ x {(0,0,0)} for any 
open set U. 

Using lemma IXfSl and proposition 3.6 of [17], we see from the above considerations that for any R > 
there exists an eo > such that the maps 

(-60, eo) X i?fl(ff|_i) X BniH^f Hi : (e, 3, tr, 0) ^ 

and 

(-eo,eo) X iHt2? x Br{HI_,) x (i/^ii)^ x BR{Hlf H^^ ■ (e, 34, /?, 3, », <^) ^ C-'" 
are analytic. Since 

C^|,=o = Am^ - 9//?^ + dai' - Kpw', C%=o = - np (5.10) 

and for —1 < (5 < the Laplacian A : ^5-2 isomorphism (see [1], proposition 2.2), we can 

use the analytic version of the implicit function theorem (see [10] theorem 15.3) to conclude, shrinking eo 
if necessary, that there exists an open neighborhood U of any point in (iJ^r2 )^ x Bji{Hg_^) x {H^Zl Y x 
Bii{Hg) and analytic maps 

(-eo,eo)x[/ — > Hi : {e, p,p,w' ' — * 

and 

(-eo,eo)xt/ — > Hi : {e, p,p,w' ,u, ' — 
such the constraints are satisfied, i.e. 

(e, p, p, ■w',u, f3, 3, "^(e, P, P, w^, 34, 3), w^^, 34, 3)) = 

for aU (e,p,p, w^,34,/3,3) e (-eo,eo) x C/. □ 



13 



6 Local existence for the Einstein-Euler system 



The combined systems (|3.19|) and (I4.12p can be written as 



b°{eU, eV)dtV = -c'diV + 6^(e, U, V)diV + /(e, U, V)V + -giV)V + K (t ^ x") (6.1) 



where 



b"{eU,eV) 



U:= (0,0,u^^O,Of , 

O 

A^{eu) 

a^(eu, ea) 

oj ' 



6u'^ := u*^ - u'^ 







f{e,U,V)V:= 



^{e, eu, w% a) 

/ef ^■(eu,Ufe) - S'^ + Ae6u'^dkp'' - 8edkP^'Su^^''\ 




V 



b{e,eu,Uk,w\a) 



I 



g{y)V:= (-^l,5ip(a),0,...,0f , 



and 



I O O 

V 



(6.2) 
(6.3) 

(6.4) 
(6.5) 
(6.6) 

(6.7) 
(6.8) 

(6.9) 



For initial data, we will use the following notation: given a function z that depends on time t, we define 



To fix a region on which the system (|6.ip is well defined, we note from (|3.20p . (|4.20p . and the invert- 
ibility of the Lorentz metric (77*-' ) that there exists a constant Kq > such that 



det(?7*^' + 4eu*^') > 1/16 , 1 + ew^ > 1/16 , 
^^(eu) > , a\eu,ew,ea) > 



(6.10) 
(6.11) 



and 



\A'^{eu)\ < 16, |a''(eu,eu',ea)| < 16 (6.12) 

for all |eu| < 2X0, |ew'| < 2Kq, \£a\ < 2Kq. The choice of the bounds 1/16 and 16 is somewhat arbitrary 
and they can be replaced by any number of the form 1/M and M for any M > 1 without changing any 
of the arguments presented in the following sections. However, since we are interested in the limit e \ 0, 
we lose nothing by assuming ill = 16. 

Proposition 6.1. Suppose -1 < (5 < 0, /c > 3 + s, a, G H^,, i^'' e ^i' € Hl_., (3^ e 

O O 

C^{[—T^T]^H^_-,). Let vJ'^ , dtvJ'J and be the initial data constructed in proposition which, by 

00 o 

choosing < 1 small enough, satisfies 

\ew'\,\ea\,\ev^^\<Ko for all e e {Q,e^]. 



o o o 



Then 
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(i) for each e G (0, eo]j there exists Ti{e),T2{e) > and a unique solution 

V,€^C'{{-T,(e),T2{e)),Hti) 
to the system (j6.ip with initial data 

Q o o o o 

(ii) The identities 

^tu^' = ^ , and u'l^ = ^jii:^' 
hold where by definition u*^ = e^^ul^ , u*-' = + <^u*-', and = eUj. 

o o o 

(Hi) The triple {ul^ ,wl,a^} determines, via the formulas (|1.12|) . (|3.4|) . (j3.5p . and (j4.1ip . a solution to 
the full Einstein- Euler system (jl.ip - (|1.2p t/iat satisfies the constraints 

edtut^ + divJJ = anrf w'ui = — . 

(^wj For some constant C > independent of e, the initial data Vg satisfies the estimate 

o 

\\Ve-Vo\\H>, <C\\V,-Vo\\h^ <Ce 

O O 5 — l,e o O 6 — 1 

while 

for all e S (0, eo]. 

(v) If s^Po<t<T2(e) ll^e(^)ll < and for all {x,t) €E x [0,T2(e)) |e5Ue(x,t)| < i^o, |ew*(a;,t)| < 
2Kq, and |eQ!e(a;,t)| < 2Kq, then there exists aT<,> T2(e) such that the solution can be continued 
to the interval (— ri(e),T*). 

Proof, (i) Follows directly from theorem IB. 5| proposition IB . 6| and corollary IB. 7| where we use the initial 
data from proposition 15. II 

(ii) This follows from standard arguments on reductions of 2"'^ order hyperbolic equations to l*** order 
symmetric hyperbolic systems. See [39], section 16.3 for details. 

(iii) By part (ii), the triplet {u^-' , w*, cte} satisfies the reduced Einstein equations (|3.17p and the fluid 
equations (|4.12p . By construction, {ul^t=o,wl\t=o,ct^\t=Q} satisfies the constraints N\t=o — 0, W\t=o = 
0, and (CJ^* — T'*')|t=o = 0. The reduced Einstein equations p.l7p can be written in terms of the Einstein 
density Q^^ as 

g'^ - -^'dkU^ + IdkUH'^^ = T"' ■ (6.13) 
Using (g4i _ T''*)|t^o = 0, we see that 

{--Q^^dkU^ + 24^(^0^)'=) = . 

A straightforward calculation then shows that this implies that dtH-'\t=o = 0. As discussed in sections |4] 
(see (|4.24l) ). M satisfies a linear symmetric hyperbolic system and hence by uniqueness, it follows that 
A/" = for all {x^ ,t) e x (-Ti(e), r2(e)- Thus {wl.tte} determine a solution, via the formulas (|4.1ip . 
to the Euler equation which are equivalent to ViT^^ = 0. So taking the divergence of (|6.13p while using 
'SJ .{f^i = WiG'^^ = shows that W satisfies an equation of the form 
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where the Q^^ are analytic in g and dkQ- Clearly, this is a linear, 2"*^ order hyperbolic equation for W . 
Since W\t=o = dtW\t=o = 0, we must have W = for aU (x^, i) e x (-ri(e), T2{e)). 

(iv) We know from proposition 15 . II that the map (0, eg] 3 e — > V^e G ^l-i e analytic which implies the 
estimate — Fo||_f/fc < C'e for some fixed constant C > 0. So then 

\\V,-V^\\h^^ < \\V,-V^\\hi <Ce 

O O 6 — I.e. o O — 1 

by lemma [A. Ill Since {u€,^^;%ae} solves the reduced Einstein equations (|3.17[) . we have that 

where the /^'^ are analytic and quadratic in duu^ and 9fcUe while S^^ are also analytic and linear in 
and 'w\. Evaluating this equation at t = 0, and using the following facts from proposition 15. II 

e^'lluf-'ll^.+i + lluril^j^. + WdtU^'WH-^ ^ + \\aA\H^^ ^ + ^ < C, (6.14) 

O S O 5 O o — 1 O o — 1 O o — 1 

we find upon solving for dtdiU^^ that 

\\dtd4U^^{0)\\j^^-i <C V e e (0,eo] (6.15) 

by the calculus inequalities of appendix [X] But from part (iii), we get that 9411^^ + djul"^ — and hence 
differentiating this with respect to t and evaluating at t = yields 

||at54£if (0)11^.-1 = Wdidtui^^.-i <C V e G (0, eo] . (6.16) 
From the estimates (|6.14p . the fluid equations (|4.12p and similar arguments as above show that 

WdtaMlln^-^ ^C+\\dtwl{0)\\^.-i <C VeG(0,eo]- (6.17) 

Estimates ((04l) - ((6Tf|) and lemma ETD then imply that \\dtVe{0)\\f^k-i < \\dtVe{0)\\f^k-i < C for all e 
G (0,eo]. 

(v) This is just a statement of the continuation principle of theorem IB. 61 □ 

7 The Newtonian limit 

Let {Ve, < e < eg} be the sequence of solutions from theorem 16.11 where we will always assume that 

— 1<5< —1/2 and suppa C Br for some R> 0. 

o 

If we let Tm{e) denote the maximal time of existence for the solution V^, then 

s+l s+1 

K G fl C\[0,T„M),Htl) C fl C\[0,T„M),HtU. (7.1) 

1=0 £=0 

So G rifco ^^([O) T'm(e)), ^5^1 ) and hence proposition 3.6 of [17] and lemma fATSl implv that 

Using proposition 2.2 of [1], we can solve the equation 

A$, = (7.2) 
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to find 

s+l 

To obtain the Newtonian limit, we use $<: to take care of tfie singular term e^^g{Ve)Ve in (|6.1[) by 
introducing the new variable 

:= (ut , <5u:^ «e, wl) := xCj^^ - 5\5idj^, . (7.3) 

Observe that 
where 

d*, := (0,5^5iaj$„ 0,0,0). 

Noting that 

b%eU,,eV,) = b°ieU,,eW,) and {e,U„V,) ^ b\e,U,,W,), (7.4) 
We satisfies the equation 

6"(e?7e, eiye)9tM/e = ^c'diW, + 6^(e, f/^, W,)diW, + /(e, C/„ M^, + d^,)W, + H, (7.5) 

where 

He := /le - 6°(eJ7e, eW^O^td*. + 6^(e, C/^, We)did<^, + /(e, t/^, M^^ + d$e)d$.. 
By construction the initial data is bounded in H^_i as e \ 0. Therefore by lemma lA.lll there 

o 

exists a constant Ki such that 

\\W,\t=o\\H^_^^^ < Ki for aU e E (0,eo]. (7.6) 

Also by definition of W^^ and lemma IA.7[ 

max{\\Su,\\L^,\\a,\\L^,\\wl\\L^} < \\W,\\ci < Csob 1 1 VF, 1 1 ^ ^ (7.7) 

where Csob is the constant from lemma IA.7I that is e independent. Shrinking cq if necessary, we can 
always assume that 

2eoCsohKi < Ko . (7.8) 

Define 

Te minjsuplr > 0| sup ||M^e(t)||^fc < 2Xi and sup HKllfffc < ooj , 1 |. (7.9) 

0<t<r 0<i<r 

From the continuation principle in theorem 16. 1[ it is clear that satisfies 

< r, < T,„(e). 

7.1 Energy estimates 

We will now use energy estimates on the Hg_^ ^ spaces to show that is bounded below by a constant 
independent of e. The strategy we use is that of [5, 19] adapted to the Hg^ spaces. All of the results 
below will be derived under the assumption that the 1-parameter family K of solutions has the additional 
regularity 

y,e nC^([0,r,],i7,^+i-0. 

^=0 
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It is then not difficult to use solution of this type to approximate solutions of the regularity type (|7.ip 
and thereby show that all of the following results also hold for solutions with the regularity (|7.ip . Since 
these sort of approximation arguments are standard, we will leave the details to the interested reader. 

The next lemma contains the basic energy estimate which is the key to deriving estimates independent 
of e. We note that this type of estimate has been derived previously for the standard Sobolev spaces 
in [5,19]. It also makes clear why we need to introduce the variables and $e to put the Einstein-Euler 
equations into the form (|7.5|) . 

Lemma 7.1. Suppose eq > 0, a" e C^{[Q,t],W^'°°), € C°{[Q,t],W^'°°), g e C"{[0,t], Ll J, and that 
w e C^([0, t], i?^ is a solution to the linear equation 

a^dtw — a^djw + g . 
Then there exists a constant C > independent o/ e G [0, eo] such that 

^(u;|aMLi., < C[{\\dWaU^+e\\dU^)\\w\\l.^^ + \\ghljwhlj 
where diva = dfO^ + dia^ and a — (a^, a^, a"^). 

Proof. Let a = a~^^~^. Then ||(7~^9jfT||ioo < eC for some constant C > that is independent of e. Using 
this, the proof follows by a standard integration by parts argument as in the proof of lemma lK4l □ 

To continue, we estimate, in terms of Ki, how much the support of can change as e \ 0. 
Lemma 7.2. 

SUpp Off {t) C Br+32Ki 

for all (t,e) G [0,t,] x (0,eo]. 

Proof. Letting , Y and Y be as in section U (see (|4.26l) ). we define 

X'At) := X\ewt{t),wHt)) 

and 

:= Y{€wt{t),eaSWdtwt{t)+diwi{t))+Y{€{n,+Sn{t)),ew\t),mk{t),ew 

Using (HHni), (EH), (EH), and d?!]), we obtain the bound 

||X/(i)||ioo < 32ifi V (t,e) e [0,r,] x (0,eo]- (7.10) 

From lemmas ^Kj^ and (IXTOll . and ([731), it follows that X^ e C°{[0,t,],CI) and F/ e C°i[0,n],C^). 
Therefore the vector field X^ can be integrated to get a flow ipl{t,x) that is well defined for all 
{t,x) G [0,n] X R^. For each x G K^ define af(i) := a,{t,ip,{t,x)). Then dt^l{t,x) = Xl{t,ip,{t,x)) 
together with the evolution equation (|4.25p implies that 

ja^,{t) + Y{t,^,{t,x))a:{t) = Q. 

By assumption suppao C Bj^ and hence 0;^ (0) — a{x) — for x G Ej^ M.^ \ Bji. Therefore 

o 

a^{t,ij^{t,x)) = allxeEji (7.11) 
by the uniqueness of solutions to ODEs. But 

Jo Jo 

by (TTTU)) and < < 1. From this, (frTT|) . and the fact that for each t the map M'"^ 9 x tp^{t, x) G 
defines a C-^ diffeomorphism, it follows that suppac(i) C -B7?,+32Ki for all (i, e) G [0,Te] x (0,eo]. □ 
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Next, we estimate ||<I>e||^fc+2 in terms of HWcH^fc ^ . 
Lemma 7.3. Let R = R + 32Ki and 



Then there exists a constant C > such that 

d d — 1 ,e 

for all {t,e) e [0, r,] x (0,eo]- 

Proof. By lemma [7^ the suppQ;e(i) C Br+^2Kx for all (t, e) e [0, r^] x (0, eo]- Letting R = R + 32Ki, it 
follows directly from the definition of the weighted norms that 

\\uu.<\\uh.^<{i+Rr'-'/'\\uu. 

for all functions u whose support is contained in Bj^ and for any e e (0, 1] and —rj — 3/2 > 0. Therefore 

\\Pe\\H^_^ < CCi\\pe\\Hl_^^^ 

where C > is a constant independent of e and 

Ci = (1 + R)-'^^-^^-^/^,Ji + ii + Ry 



\2k 



Since A : ^ is an isomorphism and A<i>e = p^, we have j:^fc+2 < C||pe||jLffc and hence, 

by lemma EH (see also (|1.12l) and (I7.3p ) and the above estimate that 

< CCillp.ll^. ^ < CCilla.ll". < CC,\m\H>; , • 

O o — d — l,e o — l,e 

□ 

We note that for the remainder of this section, all of the constants appearing in the estimates may 
depend on the fixed constant Ki . We will often use C to denote constants that depend on Ki and that 
may change from line to line. 

Let = D'^W, {\a\ > 0), 6° = 6°(e{7„ eW^), = b^{e, C/,, 14) and /, = /(e, U„ + The 
evolution equation (|7.5p implies that 



dtW, = Qc^ + bi^ diW, + + {b^r'H, . (7.12) 

Differentiating this equation yields 

b%W^ = ^c'diW^ + bidiW^ + q°' \a\ > (7.13) 

where 

= b'^[D", {b"J-\e-'c' + bi)]diW, + b'iD"{(bl)-^f,) + b"^D"{{b'J-'H,) . (7.14) 
From lemma lA.lll we know, since— 1 < S < —1/2, that ||eUj||„fc+i < el''+-^/^l ||ue|| . Since ||u£||„fc+i 

O 5.E o 6 o S 

is uniformly bounded in e, we get, by lemmas [A. 71 and I A. 1 11 that 

b S.e 

for some constant C > independent of e. So 

Wbimw-.^ <C V(t,e) G [0,T,] X (0,eo] (7.16) 
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by (HH), (1721), ([73 and ([71^ . Also, note that 

+ pd$jLoo < < C and lia^d^Jicc < 

by (|A.3|) . (|A.24p and lemmas FTSl and lA. 71 The evolution equation ((7712)) then implies that 

\\^th%L^^\\eDb\eU,,eW,)■^tW,\\L^<C{l + \\^td^,\\H.^^). (7.17) 

Together (|7.16p and (|7.17p establish the existence of a constant C > such that 

lldivfeeWIU- < + V(t,e) e [0,T,] X (0,eo]- (7.18) 

Differentiating (6°)"^ yields 

d.Mr^ = -e(6:^)-i(™0(eC/„eT^,) • (a^f/,, cJjW^,)) (6°)-i . 

This along with (|7.15p . (|7.16p . (|A.3[) . (jA.24[) . and lemmas [A . 71 and [A . 91 can be used to control the singular 
term in (|7.14p and results in the following estimate (see also appendix IB. 2^ 



lk"(i)IL2 <p„(||w^,(t)||„. ||$,(t)||^.^.,lia,$,(t)ll .+0 vte[o,T,] (7.19) 

o — i- — |a|.c o — i,c d 

where Paivi, 2/2,2/3) is a polynomial that is independent of e and satisfies P(0) — 0. Note that in deriving 
this results, we have used the estimate 



+ <c||$,||^.+2 and <q|at$,||^.+i (7.20) 

for some C independent of e which follows from (jA.3p . (jA.24p . and lemma lA. Ill 
Define 

\a\<k 

Then 

^\\WM\hI_,^^ < |||W^e(i)|||M-i,e < M\WM\h^_,^^ Vt e [0,T,] (7.21) 
by (|CT1) and (pT^ . Lemma O combined with the estimates (TTTC)) . ([7t5)) . and (fTTO)) implies that 

or equivalently 

^|||TV.(t)||U.5-i,. <P(|||M^.(i)IIU.5-i,.,ll<fe(i)llff''+^,l|5t$.(i)||^^+0 Vie [0,T,] (7.22) 

for a e independent polynomial P(j/i, ?/2, J/s) satisfying P(0) = 0. By lemma [731 ||$£||^fc+2 can be bounded 
by a polynomial of ^ that is independent of e and vanishes for || Wejj^fc ^ =0. The differential 

inequality [722] shows that if we can do the same for ||5f$e||^fc+i then we get an estimate for |||T44(/;)|||fc_5_i^(: 
independent of e. 

Lemma 7.4. There exists a polynomial P{y) with coefficients independent of e such that P{0) = and 

\\dt<i>M\H^+^ < Pi\\wM\Hi_J 

for all {t,e) £ [0, r,] x (0,eo]. 
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Proof. By (I4.12p . We :— {a^,'wl)'^ satisfies an equation of the form 

and so 
Thus 

d — 1,6 l,e d— 2,e (),e o — l,e U,e d — l,e 

by lemma [Ol Also by (fTTS)) . (|A.3p . (|A.24p . and lemmas [Ol and I A. 91 we have that 

l.c o — J-,c — o — J-,c 

||(a4)-i6i||^.-i <P(||W^e||^^._^ J, and || (a^)-!^^!!^,^, < P(|| J 

for some polynomial P{y) that is independent of e. The above two inequalities along with (|7.20p and 
lemma 17.31 show that 

for a polynomial P{y) independent of e and satisfying P(0) = 0. Using lemma [X?8l the above estimate 
implies that 

\\dtp4\H^:i^<Pi\m\Hi_J 

where as above P{y) is a polynomial that is independent of e. Since Adt^e = dtPe, the same arguments 
used in the proof of lemma 17.31 can be used to conclude 

□ 

Lemmas 17.31 and 17.41 combined with the estimate ()7.22p yield 

^|||T^eW||U,5-l,. < Pmemk.S-lMWemk,5-l,e G [0, T,] (7.23) 

for a polynomial P{y) that is independent of e and whose coefficients depend only on Ki. By Gronwall's 
inequality there exists a time T* € (0,1), independent of e, such that if y{t) > is and satisfies 
dy/dt < P{y)y, then y{t) < e^^*y{0) where K3 is a constant that depends on Ki. Therefore 

\lW,{t)\lk,s-i,e < e^^*|||P^,(0)|||fc,5_i,, for all (t, e) e [0, min{T*, rj] x (0,eo]. (7.24) 

Shrinking T* if necessary, we conclude that 

\lWemk,5-i,e <^Ki for all (t,e) e [0,minr*,re] X (0,eo]. (7.25) 

Note also that 

\\VemH^,_,,, < C for all (t,e) G [0, min{T*, rj] x (0, eq] (7.26) 

by I7.20i 17.211 and lemma 17.31 Therefore by the definition of , we must have < T* < for all 
< e < eo. 

Differentiating (I7.12p with respect to t, shows that We '■— dtWe and d$e := dtd^e satisfy the equation 
6°(eC/„ eWe)dtW, = ^c' diW.+b' (e, U,, We)diWe + /i(e, U,, W^,, W„ d*,, d$,)M^e 
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for analytic functions /i, /2 with /2 linear in the last 3 variables. This equation has the same structure 
(|7.5p and is not difficult to show that the arguments used to derive (|7.24p can also be used to obtain the 
estimate 

under the assumption that ||Vl4(0)|j ji^t-i is bounded as e \ 0. But this is clear from proposition 16 . 1 1 and 
lemma [7751 and so the estimate holds. We have proved the following proposition. 

Proposition 7.5. For eo > small enough, there exists a T* > independent o/ e € (0, Eq] such that 
the one parameter family of solutions exist, for all e G (0,eo], on a common time interval [0,T*]. 
Moreover, there exists constants C > 0, R > such that 

max{\\Su,\\L^,\\a,\\L^,\\vul\\L^}<^ \\VM\m , <C, \\dtK{t)\\ ^.-^ < C , 

<c, ||ata>,(t)||^.+i <c, 

andsuppa^{t) C Bj^ for all (e,t) G (0, eo] x [0,T*]. 

7.2 Properties of the limit equations 

To fully understand the limit equations of section [7.31 we first need to consider the following system 

dta ^ -w'dia ~ ^diw^ (7.28) 
2n 

dtw-' = - —d^a - w'diw^ - d^^ (7.29) 
2n 

A$ = /5 (7.30) 

with initial data 

d(0) = a and w\0) = (7.31) 

where a and are as defined in proposition 16.11 This system is precisely the Poisson-Euler equation 

written using the Makino variable p = {4Kn{n+i))-" O'^"- Indeed, a straightforward calculation shows that 
(p, w^) satisfy the Poisson-Euler equations of Newtonian gravity 

dtp + di{pw') = 0, (7.32) 

p{dtw-^ + w^diw'') = -{pd-'^ + d-^p) , (7.33) 
A<l = p, (7.34) 

where p = X/5("+i)/". 

Proposition 7.6. There exists a T > and a solution 

a, e C°([0, T],Ht,) n C\[0, T],Hlzl) , 
$ e C°{[0,T],H^+^)nC^{[0,T],H^+^), e C'>{[0,T],H^+l) 

to the initial value problem (|7.28p - (17.3111 where a{t) has compact support for all t G [0, T]. Moreover 

(i) this solution is unique in the class 

a,w €C°{[Q,T],H'')C^C^{«' x [0,r]) ^ C^{[0,T],H''^+^)r\C^{W x [o,r]) 

where a{t) has compact support for all t G [0,T] , and 
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(ii) the solution also satisfies 

Proof. Writing the system ((7?28)) - (fr30)) as 

I wJ) - \-i-^S" ~w' J \w') {AKn{n + 1))" {d'i^-^a'^) 



we see that this system is symmetric hyperbohc with a non-local source term. Since A : H^_2 is 

an isomorphism, it is not difhcult to adapt the approximation scheme and energy estimates of appendices 
IB.ll and IB. 21 to this system. Then as in appendix IB.3| this is enough to produce an existence theorem. 
Consequently, there exists a T > and a solution 

a,w' e C"{[0,T],Hti)nCH[0,T],Htl)- (7.35) 

Therefore 

p e C°{[0,T],Hl,)nCH[0,T],Ht^) (7.36) 

and hence $ = A-^p e C°{[0,T], H^+^) n C\[0,T], H^+^). 

Differentiating (|7.34|1 with respect to t and using (17.321) yields 

Adt^ = -di{pw^). (7.37) 

But, (17:3511 implies that pw' G C°{[0,T], HI_^) and hence A-^pw^) € C°i[0,T], H^+'^). Taking the 
divergence then gives di{A-^pw^) € C°([0, T],Hg+l). However, (TZlTTl) imphes that 5t$ = -A-^di{pw^) 
= -di{A-^{pw^)) and so dt^ S C°([0, T], iJ^^+i^). 

The statement about compact support follows from the symmetric hyperbolic equation satisfied by a 
and the property of finite propagation speed. Uniqueness follows from a slight modification of standard 
arguments, see [39] proposition 1.3, section 16.1. □ 

7.3 Convergence as e \ 

In this section, we identify the limit of the relativistic solutions as e \ 0. To accomplish this, we adapt 
the arguments of [37], section III. Define 



5/ f M 



V 

b 



A'{6v^^) 


w' 





j{V)V:= (-5^^0,il^^6)^, 



and 



h ■= {r]'^di(3^ - 2diP^\^^', 0, . . . , 0)^ 
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Theorem 7.7. For any r > 0, $e and V, converge in C°([0, T%H^+'^-'') and C0([0, T%H'^-J) as e \ 
io $ G CHK^ X [0,r*])nC°([0,r*],i7^*^+^) a«d the unique solutionV £ C\R^ x [0, T*]) n C°([0, T*], i?*^) 
o/ </ie system 

P{dtV - VdiV - /(V^)V^ - ^) = , 
c'di{V -d^) = 0, 
1^(0) -1^,(0), 

where P is t/ie projection onto the LP' orthogonal complement of { c^djW = 0\W € H^}- Moreover, 
(i) there exists a R> such that swppa{t) C Bj^ for all t G [0,T*], 
(li) there exists a lu e C°([0, T*], such that diLU e C°([0, T*], i?'^"^) and 

dtV - VdiV - f{y)V ~h- c'diLu = , (7.38) 
(Hi) and for 5i > —1/2, there exists auE C'^([0,T], L^^) such that 

Proof. By assumption —1 < 5 < —1/2, and so it follows directly from the definition of the weighted 
norms that for every ^ > 0, 

\Hh^ < \MhI_,^^ for aU u e H^,^^. (7.39) 

So by proposition 17.51 

and e C°([0,r*],i75^+2) n Ci([0,r*],i75^+^) are uniformly bounded for e € (0,eo]. Therefore by the 
Banach-Alaoglu theorem there exists subsequences of <I>e and Ve, which we still denote by <I>e and Ve, and 
$ e Li'°°([0,r*],If5^+^) nLip([0,T*],H5^), y G L1'°°([0,T*],H'=) n Lip([0,r*],i7'=-i) such that and 
Ve converge weakly to l> and V", respectively, as e \ 0. 

By proposition 17.51 the support of is uniformly bounded in e and hence the support of the weak 
limit a must also be bounded. From proposition [O] we have that Uj^ = dj\x]P . So by lemmas [A. 71 and 
and ifTTS)) . we find that for 5i>-l/2>5 

wn\Li^ < c\\um\L<i < cwK'hi^ < c{\\u'ih._^^ + i!u:^L. j < c(i + mH^_j 

for a constant C independent of e. It follows that ii{ converges weakly to a u*-' £ L^'°°([0, T*], for 
which dju^^ = uj. 
Now, satisfies 

b°ieU„ eV,)dtV, - -c'diiV, - d<^,) + b\e, U„ V,)diV, - /(e, C/^, V,)V, - h{eU,) = , (7.40) 
e 

and hence it follows from the boundedness of and that 

\\c'di{V, - < \Wdi{V, ~ d^^H'-l^ < Ce. 

Letting e \ yields 

c'diiV-d^) = 0. 
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Next, applying the projection P (note that — d<J>e e H^) to (|7.40p gives 

P(6°(eC/„ eV,)dtV, - h\e, U,,V,)diV, - f{e, C/„ V,)V, - h{eU,)) = 

or equivalently 

where we set 6^ = b°{eUe,eVe), hi = ie,Ue,V^)diV„ /, = /(e, f/„ V;)V;, and h, = h{eUe). Suppose 
ip G and let {u\v) = J^^ uv(Px be the standard norm. Then 

(V'|PfeS?(I - r)dtV,) = ((I - P)b",P2p\dtV,) (7.41) 

as P is a self-adjoint projection operator. Since the imbedding H^{Bfi) — > H'^~'^{Bb) (r > 0) is compact 
for any ball Br, and converge in C°([0, T%H^-J) and C"([0, T*],H^+^-'') to V and respectively, 
as e \ 0. Using this strong convergence and ([77T5| . we find that (I - P)6°PV' -> (I - P)PV' = in 
as e \ and hence {tp\¥b°{l - P)dtVe) -> by (|7\4T1) and the fact that ||5tVe||L2 is uniformly bounded 
in e. Therefore, we have established that 

P&^(I - P)dtVe — > weakly in as e \ 0. 

The remainder of the proof follows from a straightforward adaptation of the proof of Theorem 2 in [37] . □ 

From the block diagonal form of the matrix c^, it is clear that oj can be written as 

Using this, we can write the system (|7.38p as 

dt& = -w^dia - ^diw^ (7.42) 
2n 

dtw' = -^d'a - w'diw' - [5"{uf + feiitf ^) + 4uf ] (7.43) 

dtw^ = -w^diw^ - (af + Sijui^) (7.44) 

dtui = ASu^'diui + ASu"diuj + ifdiP^ ~ 2dil3^'if^^ - S'^ + d^to'/ (7.45) 

dtu/ = ASu^-^dju'i + diul' (7.46) 

dtSu'i = ui^ (7.47) 

djei = (7.48) 

d-^uj = 6l6iA^ (7.49) 

A* = p (7.50) 

with initial conditions 

ui^(o) = o, af(o) = o, af = a,0 (0:= A-ip(o)), (7.51) 

a{0)^a, w'^w^, (7.52) 

O O 

Equation (|7.48p immediately implies that 

uli^O, (7.53) 
and hence, by uniqueness and the fact that ^u*-' (0) = 0, it follows from (|7.47p that 

60 = . (7.54) 

Since uj = dju'^ , we get from (fOOl) that Au*-' 6l6lA^. But u*-' G and A$ e L'j_^ and so by 
theorem 1.2 and proposition 1.6 of [1], we find that u'-' e Hg for > 62 > Si > —1/2 > S > —1. Since 
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the Laplacian A : Hg^ — > -ff^^-^i injective for 62 < (see [1], proposition 2.2), we must have u*-' — ^4^4$ 
and hence 

u'j = SlSidA ■ (7.55) 
Substituting (|735)) - (|735)) into ^T^-l^^ yields 

dta — —w^dja — —diw^ , (7.56) 
2n 

dtw-' = ~—d'^a-w^diw-^ -d-^^, (7.57) 
2n 

A$ = p, (7.58) 

and 

dtw"^ = -w'diw^ , (7.59) 

9^^^*^' = 77*Ja//3^ - 2dil3''W^^ + 5*^' , (7.60) 

diuji'' = , (7.61) 

dtdi^ - a/c^f . (7.62) 

Since u;^(0) = 0, uniqueness of solutions to hyperbolic equations implies that 

w'^^Q. (7.63) 

Proposition [7?ni and (|7.56p - (|7.58p imply that {$, w^, ci} must satisfy 

a,w^ € C''i[0,T*], H^_^) n C\[0,T*l H^zl) (7.64) 

and 

!> e C°i[0,T*], H^+^) n C\[0,T*], H^+^) , (7.65) 

dt^ e C°{[0,T*],H^+l)nC\[0,T*],Hl_^) . (7.66) 

We then get from (|73T|) and fTTB^ that 

u;t^ =dt^eC\[0,T*],H^_,) (7.67) 

and 

w^'^ . (7.68) 

Equations (|73i)) and ((7:^5)) imply that 5*-'' can be written as S''^ = 2S'fdl^pw^. We then find from (|7rBD)) 
that 

u;^ = din'^ (7.69) 

where 

17*^' = A-^{r]''diP' - 2dip^'r]^'^' + 2S'-;Si^pw^) . (7.70) 

Note that 

O'J' e c\[o,r*],i75*''+i) 

since S//?^' £ CH[0, T*], H^^Ia ) and 5*^ G CmO,T*], H^Z^) by (TTIMl) . Therefore 

= diii'^ £ c\[o, r*], i?f_i) . (7.71) 

We collect the above results in the following proposition. 
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Proposition 7.8. The limit solution {V,^} from theorem \ 7. 7| satisfies 

5u'3 = u'l = = , 

* e c"([o, t\h\-^^) n ci([o, t\h\-^^) , 9t$ e c°([o, r*], i/^+i) n ci([o, t\b\_^) , 
^i^/' = ^l<5i9j$ e Ci([0,r*],i7,^+i) n ci([o,T*],i/|_i) , 

while {<b,a,w'} solves the equations (I7.56p - (j7.58p . Moreover, the lo from theorem \7.T\ is given by 

c. = K^c.}^o,...,o^ 

where 

cjI^ = SlSidt^ € C'{[0,T*],Hl_,) , 

7.4 Error Estimate 

To get an error estimate which measures the difference between the relativistic and Newtonian solutions, 
we adapt the arguments of [37] , section IV. Define 

:= - V + d$e - - euj and 7^ := - a . 

A simple but useful observation is that 

HWh^-i =||a,-(5||„. <\\Z,\\j,.-i and |K - .-1 <\\Z,\\h^^-i . (7.72) 

— 1, c o — J-^c — 1, e — 1, c o — l,e 

Lemma 7.9. There exists an e independent constant C > such that 

\\d<^e{t) - d^{t)\\^.-i + \\Dd<^,{t) - Dd^{t)\\^.-i < - ^{t)\\j^,+i < C\\Z,{t)\\jj,-i 

5 — l,e 5 — 2,e iS S — l,e 

\\dtd<f,it) - dtd^t)\\j^.-.^ < WdtMt) - dMt)\\Hi < c\\z,it)\\^.-^^ + Ce 

and 

\\da,\\H^-^^^^ < C\\ZM\^.-^^^+ Ce 

for all {t,e) £ [0, T*] x (0,eo]- 

Proof. Since the support of a^it) and a{t) are both bounded for all (t, e) G [0,T*] x (0, eo], there exists 
a e independent constant C > such that 

o — 2 o — l,e (5 — 2 

Also, A$e = pe, AI) = p, and A : Hg^^ ^s-i isomorphism, and therefore 

W'^'e-^H"^' <\\P^~ PWh"-' <G\\P^- PWh"-' <C\he\\H'-i <C\\Z,\\jj,-l (7.73) 

5 5 — 2 5 — l,e 5 — l,g (5 — l.e 

by proposition 17.51 and lemma lA.lOl From (|4.25p and (I7.56p . it follows that 7^ satisfies 

dae = -x'dn, - Yj, + [x' - w')dia + {y- , (7.74) 
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where and Y are given by g2ni). But X' = X'{ew^,wi) and = X^{0,w^) and hence 

\\X' -w'\\jjk-2 <C\\wl -w^\\jjk-2 <C\\Z^\\jjk-i (7.75) 

by (|7J2)l . (|A.24p . lemma [OOl and proposition [731 Next, 

where r(0) = and r(0) - l/(2n) = 0. Using (|7?T5l) . (|A.3p . (|A.24p . proposition [731 and lemmas IATtI 
- lA.lOl we can estimate each of the above terms as follows 

||(y(el4) - + djwi)\\^.-2 < \\{Y{eV,) - {e\\dtwt\\^.-. + \\wi\\^. . ) 

Zn 2,11 s-i.<L -'-'5-1,8 -'-'5-1,5 

\\^{diwi - diw')\\^.-. < C||Z,||^.-i^ , 

and 



Therefore 



We can also estimate X^ and Y as follows 

<C||y,||^._^^ <C, (7.77) 
< C(||;7e||^. + 11^.11 + IIStKllH-^ ) < (7.78) 

The estimates ([777^ . ([7T751) . ([77751) . ([7T7ll . ([7T7511 along with lemma [Ol imply via the equation ([7?7il) 
that 

119*7.11^^-2^ < C||Z,||^._-i^ + Ce . (7.79) 

Since A9t$e — 5p^ and A9t$ = dtp, the same arguments used to cstabhsh the estimate (|7.73p can be 
used in conjunction with (|7.79p to show 

\\dt^,-dtnm<C\\Z,\\^.-i +Ce. (7.80) 

o 5 — l,c 

Finally from (|7.73p . (|7.80p . and lemma [A. Ill we get the desired estimates 

-d$||^fc-i +\\Dd^^- Dd^\jjk-i < - l>||^k+i < C||Ze||j^*,-i 

5 — l,e (5 — 2,e (5 5 — l,e 

and 

\\dtd^, - dtdn^u-i^ < \\dt<^, - OMh^ < C\\Z,\\^.-,^ + Ce 

for some constant C independent of e. □ 
Lemma 7.10. There exists a constant C > such that 

+\\V,it)-Vit)\\j^.-i < Ce for all (t, e) e [0,T*] x {0, e^]. 
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Proof. From the evolution equation (j6.ip . we find that satisfies the equation 

b%Z, = -^c'diZ, + h{diZ, + F, (7.81) 

where 6° = h°{eU,,eV,), hi = b{e, Ue,Ve) and 

F, = -b%{d^ - d^,) - eb%uj + bi{did^ - did<^,) + ehidiuo 

- (6^ - I )dtV + [V - bi)diV + /(e, C/„ y^V; - /(f )f + K~h. (7.82) 

Using (frT5| . (|A.3p . (|A.24p . lemmas FTQI IA.7IIA.9[ and propositions [73] and [TJl we get the following 
estimates 

- I 11^.^ < Ce{\mHU^ + \mH-_J < Ce, (7.83) 
\\b°dtid^ - d^,)\\j^.-i < ||(&° - I )dt{d^ - d'i>,)\\^,-i + \\dt{d^ - d$.)||j:^^.-i 

< C(|16°-I ll^.-i +l)||at(dl>-(i$,)||^.-i <C||Z,||^.-i +Ce, (7.84) 

||e6°9«t^||^.-i^ < eC{lb\ - + l)||9tc^||^._-i^ < Ce, (7.85) 

||e6fazt^|| .-1 <C7e||6f||^.-i||az^|| .-1 <Ct\\b{\\^. <Ce, (7.86) 

- /iIIh^^-i <Ce. (7.87) 

To estimate the term hi ~ b, we first note that 

\ a^(eUe, ew*, eofe, , Qfe) — a^(0, 0, 0, w*, a) y 

where the map is analytic. Next, the estimate (j7.15p implies that 

l|u:^l!^^^-i , < lluj^ll < C + .-1 . (7.88) 

o ' O o.e o — 1,6 — 1, e 

From proposition 15.21 and lemma [A. Ill we see that the li'^'-^ can be estimated by 

o 

lluf II ... = lleCif < el^+i/^lllu^-'ll < Cel*+i/2|+i . (7.89) 

O o,e O o,e O o 

Also, from proposition 17.81 and lemma [7^ we obtain 

\\dtae-dta\\jjk-2 + - < ||Ze||^;c-i +e||w||^fc-i < C\\Z,\\f^k-i +Ce. 

5 — l,e 5 — ■ 1 , e iS — l.e iS — l.c 5 — l,e iS — l,e 

(7.90) 

The three estimates (17.88^ - (17.901) along with lemmas [A.gi and [A. 101 and proposition l7.5l and l7.81 show 
that 

\\A\ul')+A'{Snl')\\^.-^<C\\ul'\\H| +C||<5uf|| < + , 

O O 5-1, e (5 — l,e 

and 

||a"'^(eUe, ewl, eae, , Ofc) — a"'^(0, 0, 0, a)|| j:^.-! < C(e||u^^ ||^fc-i + — a|| j:^fo-i 

+ \\wl-w'\\^^-i < Ce + C||Z,||^.-i ). 

Therefore 

\\bi-b\\„>,-i <Ce + C\\Z,\\fj,-i 
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and hence 

\\ib~bi)diV\\H^-} <C\\b^bi\\^.-4DV\\H^-} <Ce+\\Z,\\J,.^^ . (7.91) 



5-1, ( 



Next, we notice that 

f{e, U,, V,)V, - f{V)V = -P.J', + f{V,)V, - f{V)V, + e/(e, U,,V,)V, 

where 

:=-4p,(6,5^5i77p,uP«,0,...,0f 

O 

and / and / are analytic. We obtain 

\\f{V,)V, - f{V)V,+eJ(e,U,,V,)V,\\H.-^^ < ^ + Ce (7.92) 

by the arguments used above. Also, the boundedness of the support of a.{t) imphes that 

< Ce\\p,v^JK'\\H';-^ < Ce\\p4^.-4u\\H.-^ < C6||p,|| . i ||u|| .-i < Ce. (7.93) 

o — l,e O o a O o o — l,e O a 

So then 

||/(e, C/„ V,)V, - /(l/)y 11^.-1^ < C\\Z,\\^.-i^ + Ce (7.94) 
by (|7:g^ and Combining the estimates ([7^- ([7371), ([TTM]) . ([7^ . and (fTTM)) yields 

Letting Zf — D^Z. and differentiating the equation (|7.8ip yields 

b°dtZ^ ^ -c'diZ^ + 6^9/ + < |a| < fc - 1 , 

e 

where 

Using the estimates above along with propositions l7.5l and l7.8l and the calculus inequalities from appendix 
[Xl we find 

\\dtZ,\\^.-^ <C||Z,||^.-i +Ce, 

o-l-|a|,e 5,6 (1-2, e <S-l,e 

and hence 

h'^Wn^-^ <C\\Z,\\^.-^ +Ce. 
Combining this estimate with the estimates 

Wdtb'i + dibiWL^ <C , \\bi\\L^<C, 

and lemma [7TT1 shows that 

<C(||Z,||^.-i +e)||Z,|| .-1 0<|a|<A:-l. 
Summing over a and using Gronwall's inequality, we get 

\\Z,{t)\\H.-^^ < C||Z,(0)||^._-i^ +Ce for all (i, e) e [0,r*] x (0,eo]. 
This estimate and (|7.90p then prove the proposition since ||Ze(0)||^fc-i < Ce by proposition 16. II □ 
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We are now ready to prove a precise error estimate for the difference between the relativistic and 
Newtonian solutions. 

Proposition 7.11. Suppose — 1<(5<— 1/2 and fc > 3. Then there exists a constant C > such that 
Wn^t) - SlSimW^e ^ + Wdju^t) - 6l6idm\\^.-.^ + \\v'{t) - ^^(i)ll^.-i^ 
+e~^\\v\t) - l\\„^-i^ + \\Pe{t) - p{t)\\H^zL + ll^*'"^^^) ~ dtp{t)\\H^zL - 

for all {t,e) e [0,T*] x (0,eo]- 

Proof. From the evolution equations and proposition [7?H1 we have 

and hence integrating yields 

6||u:^(i) - SlSimiLl, < ~ SlSm^._^ + WiCiis) - eu;l^is)h. ds . (7.96) 

But ^ ^ 

/ \\u^{s)-eiul\s)h. ds< f \\V,(s)-Vis)\\J,.-^ + e||a;'^^(s)|| (7.97) 

and 

by the calculus inequalities of appendix [X\ and proposition l5.ll Also, by lemma [X4l and u/^e = 9/ lie, we 
have 

||u:^- - SlSm^e ^ < Cllu^:, - SlSid<t\\^._^^^ + eWnnt) - SlSimW^. ^ + e\\unt) ~ 6l6m)\\^.^^ . (7.99) 

Recall that = {4:Kn{n + l))^"af^ and p = {4:Kn{n+l))^'^a^'^. Since ||ae||^fc ^ is bounded as e \ 0, 
we obtain 

\\p^-p\\^^-i < C||ae - <C\\Ve~V\\jfk-i (7.100) 

by lemma lA. 101 We also have that ||9fae|| is bounded as e \ 0, so the formulas 

(4An(n+l))" (4An(ri + l))" 

and the calculus inequalities of appendix |A] imply that 

\\dtPe - dtp\\ijk-2 < C{\\ae - Q!||^^'-2 + \\dta^ - dta\\^k-2 ) 

< Ci\\V, - VWh^zL + 1'^*"^ " dtaWjj^-^). (7.101) 
Finally, from the definition of Ve and V, we have 
\\diul'{t)~Sl6id^{t)\\j^,-i +\\v\t)~w'{t)\\jj,-i +e-^\\v\t)-l\\jj,-i < C|lV;-F|l^.-i . (7.102) 

5 — l,e i5 — l,e iS — l.c 5 — l,e 

The proof now follows as a direct consequence of lemma [7. 101 and (|7.96p - (|7.102p . □ 

In the above error estimate, the norm itself depends on e. We now show how to choose norms 
independent of e which are compatible with the error estimate above. First, for any r/ G M define a norm 

by 

\a\<l 
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Recalling that -1 < 5 < -1/2, fix 77 € [S, -1/2]. Then from (|A.24|) and lemma [OTl we get that 

lk||f,2,„-i<Ce"+i/2||^||^^^^^^^ \\u\\o,e.,<Ce'^+^/^\\u\\Li^ (7.103) 

for some constant C > independent of e. Combining (|7.103p with corollary 17.111 yields the following 
theorem which is our main result. 

Theorem 7.12. Suppose —1 < S < —1/2, ~S < rj < —1/2 and fc > 3. Then there exists a constant 
C > such that 

\\unt) - SlSl^t)\\o,e,n + I! " ^4^i5,$(t)|U_i,2,,-i + \\v'it) - 7i^(t)|U_i,2,,-i 

+e-'\\v\t) - l||fc-i,2,„-i + IIp.W - p(i)||fc-i,2,„-i + \\dtPe{t) - dtpmk-2,2,v-i < Ce^+'Z' 

for all {t,€) e [0,T*] x (0,eo]- 

Note that for j] = —1/2, we have 

I|m||o,6 -1/2 = IIuIIl" and ||w||f,2 -3/2 = 11^11^^ 

where is the standard Soboley norm. So the above theorem shows that the difference between the 

relativistic and Newtonian solutions is of order e with respect to the norms || • ||^6 and || • ||^fc-i. 



A Weighted calculus inequalities 



In this and the following sections C will denote a constant that may change value from line to line but 
whose exact value is not needed. 

Let y be a finite dimensional vector space with inner product (-j-) and corresponding norm | • |. For 

^roc( 



u £ Lf^^{W^, V), 1 < p < 00, (5 G R, and e e K>o, the weighted norm of u is defined by 



— S — n/p 



u\\lp if 1 < p < CXI 
if p = 00 



where fTe(x) := y 1 + — |e.i:p. The weighted Sobolev norms are then defined by 



Q|<fc 



1/p 



\u\\^k,p := < 



|a|<fc 



if 1 < P < 00 
if p = cx) 



where k G No, a — (ai, . . . , a„) G Nq is a multi-index and D°' = d"^ . . . 9^". Here 

d 



where (x^, . . . , x") are the standard Cartesian coordinates on M". 
The weighted Sobolev spaces are then defined as 



(A.l) 



(A.2) 



w; 



k,p 



{«G<„f(M",y)|||«|U.,. <oo} 
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Directly from this definition, we observe the simple but usefril inequality 

We note that W^'q are the standard Sobolev spaces and for e > 0, the W^f are equivalent to the radially 
weighted Sobolev spaces [1,7]. For p = 2, we use the alternate notation Hg^ := Wgf. The spaces Lg^ 
and Hg^ are Hilbert spaces with inner products 

{u\v)i^2 := [ {u\v)a:^'--d^x , (A.4) 

and 

{u\v}hI^ := ^ {D<^u\D»v}l2_^^^^^, (A.5) 

\ct\<k 

respectively. When e = 1, we will also use the notation W_^'^ = W^f and J?| = -fff i- 

Let Bfi he the open ball of radius R and qr and denote the annuli B2R \ Br and Bm \ Br, 
respectively. Let {<j)j}jlQ be a smooth partition of unity satisfying 

supp^oC-B2, supp^jCAaa-i (j>l), and (j)j{x) := (f>i{2^~^x){j > 1) . 
Scaling gives a one parmeter family of smooth partitions of unity 

cP'jix) := (l^jiex) (i>0) 

which satisfy 

supp0^ C B2/, , supp0^- C (j > 1) , and Sjcj)'j{x) := ^l{2^-'x){j > 1) . (A.6) 

Define a scaling operator by 

Sjuix) := u{2^-'^x) . (A.7) 
This operator satisfies the following simple, but useful identities 

Si = 1 , Sj o Sk = Sko Sj = Sk+j-i , (A. 8) 

SjcP'j=<Pl (j>l), (A.9) 

\\S,u\\L.=2'^\\uh., (A.IO) 

and 

Sj 0"= 2(i-J)l«l£)« o Sj . (A.ll) 
Lemma A.l. For 1 < p < 00, there exists a constant C > independent 0/ e > such that 

{j o,€ ' * o,e 

Proof. From the identity 

" \u\Pd''x 



Bi/f j—1 J a^j+i 



and a simple change of variables, it follows that 
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This identity and 



-Sp-n _ J 2 ''2 " if -5p - n > 



niax a,{x) , 

2:6-84/. I 1 if -5p-n < 

, . X f 1 if -Sp -n>Q 

2 if -Sp-n < 



(1 + 22J)^^ if-Sp-n>0 



a:ea^/^ J ; | ^ ^ 22^- 1) ) if -^p - n < 

mmj5,a,)(:.) -| (1 + 2^^)-^ if -<5p - n < 



show that 



and 

L2-P'U'''^\\S,u\\U^,,^) < 2"(^-^'||5,(a7^-«/^'^.)|U.(,,^^) < C2-P'<-^~'^\\S,u\\Ua,,,^ (A.14) 

for some constant C > which is independent of e > 0. Using a change of variable, the inequahty (|A.14p 
can be written as 

^2-^'^0--i)2(i-^>l^r^,(^^^.^^^ ) < 2"(^-^'||5,(ar^-"/^«)|U.(a4/.) < C2-P'^^-'h^'-^>\\u\\l^^^^^^^^ ^ . 

(A.15) 

From 

2 2 
E'^feU4/==I^v. and E'^;+'cL^^, =Ia,,,,,^ (A.16) 

fc=0 fc=0 

and (|A.10[) . we obtain 

lk|IL(B4/j < CiUluWlp + ^ , (A.17) 

fc=i 

and 



Iwir.^.^ . < C^2"(^+'^-)||5,+fe(<^j+,z.)||i, . (A.18) 



IIP 

'^"(02^ + 1/ J 

fc=0 



Combining (|XT2)) with the inequaUties ((XT3)) . (|XT5|) . (fOTl) and (fOS)) yields 

00 

Ikll^j < C{Ut>u\\lp+Y,2-P'(^-'^\\SMju)\\U (A.19) 

for some constant C > independent of e > 0. 

Since supp^Q C -B4/e and ||0o||l^ = ||</'o||l~, we get from (IA.13|) that 

UouWl. < UoWUMUb,,,^ < C\K''-^MUb,,,-, (A.20) 
for some constant C > independent of e > 0. Next, 

2-P^O-i)||5^.(0^^„)||P^ < 2-J'*'(^-i)2"(i-^)||0j7.|U.(^^^._,^j by (Em and M), 

< 2-f^(^-i)2"(i-^)||<^i|U^|l</>>|L,(ui_^,^^_ J since \\Ml^ = ll'/'jlU-- 
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So there exists a constant C > independent of e > such that 

r c(\\u\\p,f, ^ + ||u||^p, ifj = i 

Therefore 

oo 



>2 ■ (^-21) 



IP 

I LP 



< by (imi) 



where C > is a constant independent of e > 0. The proof then follows from this inequality and 
The above lemma shows that the norm 

oo 

lll^lll^. := ||<^^u||^p +^2-^^(^-i)|15,(^»|!ip 

is equivalent for 1 < p < cx), independent of e > 0, to the weighted norm ||w||^p . For p = oo, the 
appropriate norm is 

:=sup{||0^u||l»,2-^(-''-i)||0>||l- (j > 1)} 
and it is easy to see that there exists a constant C > independent of e > such that 



^MLt, < ML^,<C\\u\\Lf;^ 

The same arguments used in proving the previous lemma can be used to establish the following 
generalization. 

Lemma A. 2. For 1 < p < oo. let 



and for p = oo let 

M^..^ :=sup{||(/.;^u||M.'e.=o,2"*(^'-i)||5,(0»|lM.-.e (j > 1)}. (A.23) 
Then there exists a constant C > independent o/ e > such that 

For the remainder of this section, we will use the two equivalent norms || • H^fc.p and ||| • l^yfc.p in- 
terchangeably and refer to both using the notation || • ||y(/fc,p. From (jA.22|) . it follows that there exist a 
constant C > independent of e > such that 

||u||^fc2,p < C||u||^fci,p whenever fc2 < ki and di < 62- (A. 24) 
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Thus we have the inclusion Wg^f C Wg^f for ^2 < and 6i < 62- The representation (jA.22[) is 
particularly useful for extending estimates from the usual Sobolev spaces Wg'^ to the weighted ones 
Wgf (e > 0) as the next lemma shows. It also makes clear the philosophy behind deriving weighted 
Sobolev inequalities which is to derive global estimates from scaling and local Sobolev inequalities [1] . 
We remark that the norm ||| • |||^fc,p, as an alternate representation for the standard weighted norms 

|j • ll^^fc.p, was introduced by Maxwell in [25]. There he used the norm to define the weighted Sobolev 
spaces for non-integral k (see also [4]). Here we will only be interested in integral k. 

Lemma A. 3. Suppose > and for all u G C°°(M", V), u i—s- Fi(u) is a map that satisfies 

1 

S,F^{u) = 2-^J-^^^F^{Sju) (j>l), 
where the Fa {a = 2,3,4,5) are linear operators on V. 
(i) If there is an estimate of the form 

\\Fi{u)\\^.,,,, <C^\\F2{u)\\ 

where pi > P2, then 

\\F,{u)\\^.,.., <C\\F2{u)\\ 

for some constant C > independent of e € [0, eq] provided 5i + X > 62- 
(ii) If there exists an estimate of the form 

\\Fi{u)\\^,k„,, < Ci\\F2{u)\\w,.2.P2\\F3{u)\\wk,.,, +C2||F2( 

where — = 1 = 1 (1 < Pi < Pa < 00 a = 2, 3, 4, 5), then 

Pi P2 P3 P4 P5 

||^^iHIU?i,Pi <C{Ci\\F2{u)\\ w'=2 .P2 \\F3{u)\\^.s...+C2\m 

Si,c 52, e e 54,6 55,1: 

for some constant C > independent 0/ e G [0, eq] provided Si + X> max{(52 + S3, Si + S^}. 

Proof. We only proof part (ii) for I < pa < 00. Part (i) can be proved in a similar manner using the 
inequality 

iJ2 "P'^" - £ "P'^' > and < g < p (A.25) 

j 3 

instead of Holder's and Minkowski's inequalities. See also the proof of theorem 1.2 in [1]. 

Recall Holder's and Minkowski's inequalities which state that forl<p<g<r< 00, 1/p — l/q + l/r 
and any two sequences Oj ,bj > that the following holds 

(i:«K)"'s(i:"i)""(i:';)"' (*-^») 

j J 3 

and 

(E(",+M')"'<(i:<.?)"'+(E »?)"'■ (A.2T) 

j 3 3 
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Next, suppose j > 2. Then 



( E 



fc=-i 



J. 

k 

1 



fe=-l 



1 

fc=-l 



J. 

fe 



fe=-l 



where C > is a constant independent of e > 0. Note that in deriving this, we have used the fact that 
VF''i'°= is bounded for e G [0, eo]- From the above inequaUty, we see that 

1 1 



1 



fc=-i 
1 



2-«5(j-i) 



E ll-^5(5j+fc((?!)j-+fcM))|!H/'=5.P5 



Pi 



k=-l 



where we have used (5i + A > niax{(52 +63,64 + 65}. The above inequahty along with (|A.26p and (|A.27p 
imply 



00 



< 



1/P2 



+ C2(^2-^^^^(^-i)||F4(5,(</)>))r^,,,,,)'^''(5]2-^=^=(^-i)||^5(5,(^ 



VF''3.P3 



1/P3 



1/PE 



and hence 



Similar arguments show that 



+ C2\\F4{u)\\^P^,,.4F5{u)\\^P,.k,) . 



(A.28) 



ii</.^i^i(u)r^,,„, + ii5i(</.'ii^i(u))rv„j < c(ciiii^2(u)ii^.„.,ii^3Nii^^P3- 



+ C2\\F4{u)\\^,,..,\\F5{u)\\ (A.29) 



for some constant C > independent of e G [0, ep]. The proof now follows from the two inequalities 
(1X28)1 and (IX29ll . □ 

The next lemma is a variation of the previous one and can be proved in the same fashion. 
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Lemma A. 4. Suppose eo > and for all u G C°°(R"', V), u > Fiiu) is a map that satisfies 

<PlF^{u)^<i>lF^{[<i>l + c^\)u), 
1 

cj,)F, (u) = ( ^ 0' > 1) ' 

where 

F2 = DP2 , F3=P3, Fi = DPi , and F5 = P5 , 
and Pa (a = 2, 3, 4, 5) are linear operators on V . 

(i) If there exists an estimate of the form 

where pi > P2, then there exists a constant C > independent of e Cz [0, eo] such that 
\\Pi{u)\\y^,k^.Pi < C(||F2(M)||^fc2,P2 +e||F2u|lvi/''2,P2) 

provided 5i + X > 62- 
(ii) If there exists an estimate of the form 

\\Fi(u)\\^.„,, < Ci||F2(m)|| \\F3{u)\\w.,.P, +C2\\F2{u)\\ 

where — = 1 = 1 (1 < Pi < Pa < 00 a = 2, 3, 4, 5), then 

Pi P2 P3 Pi P5 

\\Pi{u)\\^k,,p, < C{Ci{\\F2iu)\\^k2,p., + e\\P2u\\^k2,P2)\\F3{u)\\^k.,,p-, 

+ C2(||F4(m)||^;^,!=4,p4 + e\\P4u\\^ki,Pi)\\F5{u)\\^k^,p.^) 

(54 — l,e 54, e "55 -e 

for some constant C > independent 0/ e G [0, eo] provided Si + X > max{(52 + S3, 5^ + S^}. 

Remark A. 5. By using the generalized Holder's inequality, part (ii) of lemmas \A.3\ and \A.4\ can be extended 
in the obvious fashion if there exists estimates of the form 

\\Fi{u)\\wki.Pi < C||F2(u)||^yfc2,P2 WFsi \\Fn{u)\\ 

where ^ = J^iL^ - (1 < Pi < Pi < 00), Fi is as in lemma lX^ and Fi [i > 2) are of the form Fi — Pi 
or Fi — DPi with Pi a linear operator on V . 

We will now use these two lemmas to extend various inequalities from the standard Sobolev spaces 
to the weighted ones. All of these inequalities have been derived before by various authors, see for 
example [1,4,7,8,25,30]. The new aspect here is that we show that the constants in the inequalities are 
independent of e > and hence we find inequalities that interpolate between the weighted (e > 0) and 
the standard ones (e = 0). We begin with a weighted Holder inequality. 

Lemma A. 6. Suppose eo > Q, Si — 5i + S2 and — = \ . Then there is a constant C > 

Pi P2 P3 

independent 0/ e G [0, eo] such that 



\uv\\^Pi < C||u||^P2 ||i;||^P3 

01, e S2,c 63, t 



for all u e L^^ ^ and v € L^^ ^ . 
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Proof. Follows directly from Holder's inequality and lemma IXTSl □ 

Next, we consider weighted versions of the Sobolev inequalities. 
Lemma A. 7. 

(i) For eg > and k > n/p there exists a constant C > independent of e £ [0, Eq] such that 



for all u G Wgf. Moreover u € C^^ and for e > 0, uix) = o(|a;|'') as \x\ 



oo. 



(ii) For eg > and 1 < p < n there exists a constant C > independent of e G [0, eg] such that 

< C{\\Du\\LP_^ ^ +e\\u\\LpJ 

for all u € W^'J . 

Proof, (i) The estimate < C||u|j^yfc,p for some constant C > independent of e > follows from 

the usual Sobolev inequality ||u||loo < C||u||i4/fc,p [k > n/p) and lemma [A73l Since || • ||^yfc,p for e > 
is equivalent to || • ||^fc,p, the statement u{x) — o(|x|*) as |a;| — > oo for e > follows from theorem 1.2 in [1]. 

(ii) Follows from lemma [5^41 and the Sobolev inequality ||m||^„p/(„-p) < C||-Dii||iP which holds for all 
u e W^'P where 1 <p < n. □ 

In addition to the Sobolev inequalities, we will also require weighted versions of the multiplication 
and Moser inequalities. Wc first consider the multiplication inequalities. 

Lemma A. 8. Suppose eo > 0, 1 < p < oo, fci,/c2 > fca, fca < ki + k2 — n/p, 5i + §2 < S^, and 
Vi X V2 — > V3 : (u, v) 1-^ uv is a multiplication. Then there exists a constant C > independent of 
e e [0, eo] such that 

||uw||^<=3,P < C\\u\\^^k,.p\\v\\^k2.p 

for all u G W^^-^ and v G M^^'^. 

Proof. This proof does not follow directly from lemma [A. 3|, but can be proved in a simlar fashion. To 
see this first recall the Sobolev mlutiplication inequality 

||uu||vyfc3,p < C||M||vyfci,p|lw|lvi/fc2,p (A. 30) 

which holds for 1 < p < 00, fci, k2 > k^, and ^3 < fci + fc2 — n/p. So 



o-i 



J-1 

00 



J=l fe=-l 



1/p 



00 1 



j=i k=-i 

00 

= 1 i=i 

< C||u||^;=i,p||i;||^fc2.p 
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where in deriving the third, fourth, and fifth fines we used ((X25l) . 1X30)) . and (|X26)) . respectively. □ 
Lemma A. 9. 

(i) If eq > and Si > niax{(52 + <53, (54 + (55}, then there exists a constant C > independent of e G [0, eg] 
such that 

for all u e Hl^^ D Lf^^^ and v E H^^^ D 

(ii) If eq > and 5i > niax{(52 +(53, 64 + 65}, then there exists a constant C > independent of e £ [0, eo] 
such that 

for all \a\ <k,uE H^^^ n W^;^ and v G H^-^ n i^^^,. 

(Hi) Suppose eo>0, F e C^(F,R'") ts a map that satisfies DF e Cb"^(V',M™), and 1 < |a| < k. Then 
there exists a C > independent of e £ [0, eo] such that 

\\D'^F{u)hi , , <C\\DF\\^.-4u\\'l-J{\\Du\\^.-. +e|juL. ) 
for allue Hl^nL°°. 

(iv) Suppose eo > and F E C^(V,M'"). Then there exists aC > independent of e E [0, eo] such that 
for alluE Hl^nL'^. 

Proof. Inequafities (i) — (iv) foUow directly from ()A.24p . lemmas [A . 3 1 and FA . 4[ and the following standard 
Sobolev inequalities: 

(i) \\uv\\Hk < C{\\u\\h4v\\l-- + ||u||//;c||u||loo) for aU u E H'' n and v E r\ 1°° . 

(ii) \\[D'^,u]v\\l2 < C{\\Du\\„k-i\\v\\L'^ + \\Du\\l^\\v\\hI'^-i) for aU |a| < k, u E H'' n W^'°° and 
V E H^-^ n 

(iii) Suppose F E C^{V,Mr) is a map that satisfies DF E C'^^^{V,Mr) and 1 < |a| < k. Then 

\\d"F{u)\\L2 < C\\DF\\(^k-i\\u\\^^2\\Du\\Hk-i for aU u E H''nL°°. 

(iv) Suppose F e Cb*^(F,R''"). Then \\F{u)\\h>^ < C\\F\\c^{l + \\u\\'i-^^)\\u\\h>^ for all w e -ff'^' n L°°. 

Note that we have used || • — \\ ■ □ 

In addition to the Moser inequalities, we also need to know when the map u 1— > F{u) is locally Lipschitz 
on Hg. 

Lemma A.IO. Suppose eq > 0, F E C^{V,R), F(0) = 0, 6 <0, and k < £, and k > n/2. Then for each 
R > there exists a C > independent of e E [0, eo] such that 

\\F{ui)- F{u2)\\„k^ <C\\ui-U2\\hI^ for allui,U2 E Br{hIJ. 

Proof. See the proof of lemma B.6 in [30]. □ 
We conclude this section with a lemma comparing the norms ]] • j]^? and H • H^p . 
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Lemma A. 11. 

(i) If 5 < —n/p, 1 < p < cxD, and < e < 1, then 

for all u G L^. 

(ii) If —n/p < 5, 1 <p < oo, and < e < 1, then 

INUj<NUj^<e-^-"/^'||u|U. 

for all u dz Lg. 

Proof, (i) By assumption < e < 1, and so we have eai{x) < a^ix) < ai{x) for all a; e M". By 
assumption —S — n/p > and so we get ^-^-^IVfj~^~'^l'P < < a^''""^. Therefore, directly from 

the definition of the weighted norm, we find e~'^~"^^||'u||i,p < HwUlp < ll'^llif • P^'rt (") is proved in a 
similar fashion. □ 

B Quasilinear symmetric hyperbolic systems 

In this section we establish a local existence and uniqueness theorem for a particular form of quasilinear 
symmetric hyperbolic system on the weighted Sobolev spaces In [30], we proved a local existence 
and uniqueness theorem for quasilinear parabolic systems on the spaces by adapting the approach of 
Taylor [39] (see theorem 7.2, pg 330, and proposition 7.7, pg 334) which is based on using mollifiers to 
construct a sequence of approximate solutions and then showing that the sequence converges to a true 
solution. Here, we will again follow the same approach for quasilinear symmetric hyperbolic systems and 
adapt the local existence and uniqueness theorems of Taylor (see proposition 2.1, pg 370) to work on the 
weighted Sobolev spaces. We will only provide a brief sketch of the proof since the proof is very similar 
to the one in [30] and the details can easily be filled in by the reader. Related existence results have been 
derived independently in [4] using a different method. 

The hyperbolic equations that wc will consider are of the form 

h^\u,v)dtv = b\u,v)div + f{u,v)v + h, (B.l) 
v\t=o = vo (B.2) 

where 

(i) the map u = u{t, x) is ]R'"-valucd while the maps v = v(t, x) and h = h{t, x) are M^-valued, 

(ii) b°,V,f GC^iW xR"',Mm^m) (j = l,...,n), 

(iii) b° and V {j = 1, . . . ,n) are symmetric, and 

(iv) there exists a constant a; > such that 

&°(6,6)>wlmxm forall(a,6)e]R''xffi™. (B.3) 

B.l Galerkin method 

Let j e C^(M") be any function that satisfies j > 0, j{x) = for |a;| > 1, and J^^ j{x) d^x = 1. 
Following the standard prescription, we construct from j the mollifier jn{x) := 'n~"j{x/r]) {t] > 0) and 
the smoothing operator 

Jr,{u){x) := je * u{x) = / jr,{x - y)u{y) d^y . 
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Following Taylor ( [39],Cli. 16, sect. 1 & 2), we first solve the approximating equation 

6°(w, JnVn)dtVr, = {u, J^w^)5j J,,w,, + J,,/(u, Jr,v,f)J,jV^ + Jr,h (B.4) 
Vn\t=ii = vo, (B.5) 

and latter show that the solutions converge to a solution of (jB.ip - (|B.2p as ?/ — > 0. 

Proposition B.l. Suppose Ti,T2 > 0, r] > 0, 5 < j < 0, k > n/2, vn G , u G C°([-Ti, T2], i/^), 
and h G C^{[—Ti,T2],Hg) for some T > 0. Then there exists a T» > (T* < Ti,T2) and a unique Vjj G 
C^{{—Tif,T^),H^) that solves the initial value problem (|B.4[) - (IB.5p . Moreover if sxTp^^i^rp^ < 00 

then there exists a T* G {T^,,T2] such that Vrj extends to a unique solution on {—T^,,T*). 

Proof. Fix r/ > and define 

F{t,v) := {b'^{u, JriV))^^ {Jrib-' {u, Jj^v)diJriV + J,^f{u, JriV)JriV + Jrih) . 

Then the approximating equations (|B.4p - (jB.5p can be written as the first order differential equation 
V — F{v) ; w(0) — vq on Hg. If we can show that F is continuous and is Lipshitz in a neighborhood 
of Wo in Hg, then the proof follows immediately from standard existence, uniqueness, and continuation 
theorems for ODEs on Banach spaces. 

To prove that F is locally Lipshitz, we first prove the following lemma. 

Lemma B.2. Suppose S < -f < 0, £ > n/2 and that f G C^^(R'" x R'",M,„xm). Then for each u e 
and R > there exists a constant C > such that 

\\f{u,Vi)vi - f{u,V2)v2\\Hl < C'll^'l - V2\\hI 

for all vi,V2 G Br{HI). 

Proof. Let f(0,0) = c and g{x, y) = f{x, y) - c so that 5(0, 0) = 0. Then 

/(m, vi)vi ~ f{u, V2)v2 = c{vi - V2) + {g{u, vi) - g{u, V2))vi + g{u, V2){,vi - V2) . 
Since 7 < and £ > n/2, we get from lemma IXfSl that 

\\f{u, vi)vi - f{u,V2)v2\\Hl < C(l + \\9{u,V2)\\h{,)\\vi - V2\\hI + \\vi\\hi\\9{u,vi) ~ 9{u,V2)\\h'_^ . 
By lemma lA.lOl lemmas IA.7I and IA.9[ and (|A.24p , we get from the above inequality that 

\\f{u,Vi)vi ~ fiu,V2)v2\\Hl < C{\\u\\HI!,\\vi\\HlA\v2\\Hl)hl "W2||_f/| 

where -P(2/i, j/2, 2/3) is a polynomial. This proves the lemma. □ 

Using lemma A. 7 of [30], it is not difficult to prove the following variation of the above lemma. 

Lemma B.3. Suppose 6 < -f < 0, r] > 0, £ > n/2 and that f G C^iR"" x K",M™xm). Then for each 
u G and R > there exist a constant C > such that 

\\f{u, JjjVi)DJjjVi - f{u, JjjV2)DJ,^V2\\Ht < C\\vi - V2\\hI 

for all vi,V2 G Bn,{H^). 

The proof now follows easily from the above lemmas, lemma A. 7 of [30], and the estimates of appendix 
lAj which show that for any R > the map F : ([— Ti, T2] x Bii{Hg_^) Hg_^ is continuous and moreover 
there exists a constant C > such that \\F{t,vi)~F{t,V2)\\H'' ^ W2||_f/fc forallui,'y2 G Bji{Hg). □ 
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B.2 Energy estimates 

Fix k > n/2 + 1. By proposition IB. 11 we have a sequence of solutions € {[—T{r]),T{r])], Hg) 
(0 < T(?7) < Ti,T2) to the approximating equation (|B.4|) - (|B.5|) . The goal is to derive bounds on in 
the spaces independent of rj. To do this, we use energy estimates which we now describe. 

Lemma B.4. Suppose a" £ C\[0,tIW^'°^), G C°([0, r], T^i'°°), / G C°([0,r],L2) and that w G 
C^{[0,t], L'j^) satisfies the equation 

a^dtw = Jjja^ dj JjjW + g . 
Then there exists a constant C > independent of rj > Q such that 

j^{w\a'>w)L2 < C[{1 + ||diva||i^ + ||a||L-) Hu^Hij + ||.g||LilhllLi] 

where diva = Sja" + dja^ and a — (a^, . . . , a"). 
Proof. First, we have 

— {w\a'^w)]^2 = 2(w\a'^dtw) + {w\dta^w) 1^2 — 2(w\Jria^ djJ„w) 1^2 +2{w\g)i^2 + (w\dta^w) 1^2 . 

f~j~f A A A AAA 

Letting denote the adjoint of with respect to the inner-product (jA.4p . we can write the above 
expression as 

j^{w\a''w)Li = 2{jlw\a^d,J^w)Li + 2(w;|.g>Li + ("'|5t«°^>L5 ' (^-6) 
Integration by parts shows that 

{J^wla^ dj JrfW) 1^2^ = - {dj J^wla^ JrfW) 1^2^ - {J^^wKdja^ + p~^djp)Jrjw) 1^2^ (B.7) 

where p = a^'^^^^\ Since \\p^^djp\\L=° < 00, together lemmas B.7 and B.8 of [30] and (|B.7p imply that 
{jlw\a^djJ^w)Li < ^{d,J^w\a'Jr,w)Li + C{\ + \\d,a'\\L^ + ||a||L-)||u;||i2 . (B.8) 
Again integrating by parts and using lemma B.8 of [30], we find that 

- {djJ^w\a^ J^w)l2^ < C(l + ll^jfliLoo + ||a||L-)||w||^2 . (B.9) 

The proof now follows from the Cauchy-Schwartz inequality and equations (|B.6|) . (IB.8|) . and (IB.9p . □ 

Let v'j^ — D"v,^, 6" — b'^{u, J,jV,j), ^ = 6^ (u, J,,w,,)and = /(u, J,,t;^) J,,t;^. The evolution equation 
(|B.4|) implies that 

dtvn = ib"^)-'Jnbi,djJ^Vn + {hl)-^f^ + {b^-^h . (B.IO) 
Differentiating this equation yields 

b^dtv"^ = J,V^djJ„v'^ + ff" (B.U) 

where 

ff" = 6^[i^", (6°)-i J,6^,]9, J,i>, + &",i?"((6^)-i J,/„) + b"^D"{{b"^)-'Jr,h) . (B.12) 

To simplify the following estimates, we will assume that &'(0,0) = 0. It is not difficult to treat the 
case where V {0, 0) ^ 0. Recalling that (5 < 7 < and k > n/2 + 1, we get from the calculus inequalities 
of appendix [Al and lemma A. 7 of [30] the following estimate 

l|6°[i?",(&°)-V,6^;]a,v,||i._^^^ < \\bl\\L^\\[D'^,{b^,)-'j^bi,]d,j^v,\\L2_^^^ 

< C{\\{b"y^Jvb^v\\Hl;\\9]JvVv\\Lr^, + \\Ky^Jr,m\w^.^\\djJr,vJ^k-l) 

< C[{1 + (||u||loo + \\v,,\\l--)''~^){\\u\\h^ + \\v,^\\H^)\\vrj\\wi.'^ + (1 + + \K\\w^:^)\\vv\\hi] 
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where C is independent of 77. By the Sobolev inequaHty flemma lA.Tp we have 
and hence 

for a 77 independent polynomial P{yi,y2) ■ The other terms in can be estimated in a similar fashion 
to get 

\\9"\\^_^^^ < P{\\u\U^,H\\H^A\h\\Hi) (B.13) 

where as above P(?/i, 2/2, 2/3) is an 77 independent polynomial . It can also be shown using the calculus 
inequalities and (|B.10[) that 

||div6||i^ < P{\\u\\Hi;A\vjH;^A\h\\H0 . (B.14) 

Finally, we note that 

||6||loo < C. (B.15) 

Next, if we define 

\a\<k 

then by (|B.3|) and (jB.lSp there exists a constant C > independent of 77 such that 

C-A\v,\\H^<l\vJk,s<C\\v,\\Hl- (B.16) 

Since supg^^^y < 00 and supQ<f<;y] \\h{t)\\Hk < 00, lemma (|R4| and ((RT3)) . ((RM)) . (|RT5)) . 

and (|B.16P imply that 

^llk.||||,<C(|||7;,|||,,,)|||7;,|||fc,, (B.17) 

or equivalently 

^||b.||lM<^(llk.lllM) 

for a polynomial P{y) with positive coefficients that are independent of > 0. By Gronwall's inequality, 
(|B.17[) . and proposition IB.li this implies that there exists constants T^,K > 0, both independent of 
?7 > 0, such that r(77) > and 

sup <K. (B.18) 

0<t<T, 

Using the time reversed version of the equation (i.e. sending t ^—t —t) we also get, shrinking T» if necessary, 
that 

sup \\Vr,{t)\\H>' < K . (B.19) 
-T.<t<0 * 

Finally, from (jB.lOp . (jB.lSp . (|B.19p . lemma A. 7 of [30], lcmmas r01 and lA.9i and (|A.24|) . we see, increasing 
if necessary, that 

sup \\dtvr,{t)\\jj,-i <K. (B.20) 

-T,<t<T, * 

B.3 Local existence and uniqueness 

To get local existence following the approach of Taylor (see theorem 1.2, pg 362 in [39]), we let 77 \ 
and use the bounds (|B.18P - (|B.20p obtained from the energy estimates to extract a weakly convergent 
subsequence of which has a limit that solves the initial value problem (|B.ip - (|B.2p . Since the proof is 
very similar to that of theorem B.2 in [30], we omit the details. 
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Proposition B.5. Suppose Ti, > 0, 5 < 7 < 0, fc > n/2 + 1, vq e H^, u £ C°([-Ti, T2], iJ^) 
and h e C"{[-Ti,T2],Hl). Then there exists a > (T, < Ti,T2) and a v e i°°((-r», T»), iff ) n 
Lip((-T,,T»),i/^~^) that solves the initial value problem (|B.1|) - (IB.2|) . 

Using the estimates of appendix B of [30] and of appendix 1X1 and I B . 2 1 of this paper, it is not difficult 
to adapt the proofs of propositions 1.3-1.5, pgs. 364-365, in [39] to get the following theorem. 

Theorem B.6. The solution v from vrovosition[K5\is unique in L°°((-Ti, T2), iJj^J n Lip((-ri, T2), H^^^^) 
and satisfies the additional regularity 

Moreover, if T^ < T2 and supQ<y<2^^ < co, then there exists a T* £ (T*,T2] such that the 

solution can be extended to a solution of (IB.l[) - (|B.2p on (— T,,r*). 

For linear systems, the energy estimate (see lemma 17.11 with e = 1) ensures, via the continuation 
principle of the above theorem, that the solutions can be continued as long as the functions u{t) and h{t) 
are defined. 

Proposition B.7. Suppose Ti, > 0, 5 < 7 < 0, fc > n/2 + 1, vq e , u G C°([-ri, T2], i/.^) and 
h G C^{[—Ti,T2\,H^). Then the initial value problem 

b°{u)dtv = V{u)d,v + f{u)v + h, (B.21) 
v\t=o = vo (B.22) 

has a solution 

veC"i[-TuT2],H^)nC\[~TuT2],Ht'))- 
that IS unique m L°^{{~Ti,T2), H^^J Lip((-ri, T2), ili^;^). 

Let [n/2] denote the largest integer with [n/2] < n/2 and fco — ["-/S] + 2. Then differentiating the 
solution from theorem IB. 61 with respect to t, and using proposition IB . 71 yields the following result. 

Corollary B.8. Suppose k = ko + s, u £ flLo ^^([-^1- ^2], i?^"^) and h G flLo '^^(h^i' ^2], iJ^^O- 
Then the solution from \B.5\ satisfies the additional regularity 
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